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geometrical  spreading.  Limited  data  suggest  strain  and  frequency-dependent  attenuation,  similar 
to  the  results  from  strong-motion  data  from  Salmon  and  Sterling. 

M 

c.  EFFECT  OF  SCALED  DEPTH  ON  REGIONAL  PHASES 

Spectral  characteristics  of  the  regional  phases  Pn,  Pg,  and  Lg  are  studied  for  several  NTS 
explosions  covering  a  large  range  of  yields  and  scaled  depths.  A  new  parameter,  Effective 
Scaled  Depth  (ESD),  defined  as  the  ratio  of  shot  depth  to  the  estimated  elastic  radius,  is  intro¬ 
duced;  it  is  a  measure  of  scaled  depth  that  has  been  corrected  for  the  variation  in  overburden 
pressure  with  shot  depth.  Data  from  a  common  recording  station  indicate  the  mean  Pn/Lg  and 
Pg/Lg  spectral  slopes  to  be  most  dependent  on  ESD,  somewhat  less  on  scaled  depth,  and  the 
least  on  shot  depth.  The  spectral  shapes  of  Pn  and  Pg  are  significantly  more  dependent  on  ESD 
and  scaled  depth  than  on  shot  depth.  The  enrichment  of  higher  frequencies  due  to  increased 
ESD  (or  scaled  depth)  is  the  largest  in  the  spectra  of  Pn,  somewhat  less  for  Pg,  and  the  least 
for  Lg.  The  low-frequency  asymptote  and  the  area  under  the  initial  pulse  of  deconvolved  Pn 
also  show  much  greater  dependence  on  scaled  depth  and  ESD  than  on  shot  depth. 

A  comparison  of  inter-event  spectral  ratios  of  Pn,  Pg,  and  Lg  for  shots  with  nearly  equal 
yields  but  different  ESD  or  scaled  depths  shows  significant  differences.  Similar  variations  are 
observed  if  PSRV  response  spectra  based  on  Pn  and  the  entire  seismogram  are  compared. 

These  results  imply  that  scaling  of  the  nuclear  seismic-source  function  may  vary  significantly 
from  one  regional  phase  to  another.  Therefore,  scaling  relationships  based  on  the  entire  seismo¬ 
gram  (e.g.  Mueller  and  Murphy,  1971),  generally  dominated  by  the  Lg  wavetrain,  may  not  be 
appropriate  for  Pn. 

D.  DECON V OLUTION  AND  SPECTRAL  RATIOS  FROM  SOVIET  SALT  SHOTS  f 

The  maximum-likelihood  multichannel  deconvolution  method  of  Shumway  and  Der 
(1985)  is  applied  to  teleseismic  P  waves  from  ten  salt  shots  in  the  Astrakhan  region  of  the 
USSR.  The  data  used  came  from  the  GBA  array  and  the  four  RSTN  stations,  RSON,  RSNT, 

RSNY,  and  RSSD.  The  deconvolved  source  functions  for  each  shot  show  clear  P  and  pP 
arrivals  separated  by  about  1  sec,  suggesting  large  shot  depths.  There  are  also  several  arrivals 
which  consistently  show  up  between  the  P  and  pP,  similar  to  those  observed  for  Salmon. 

These  earlier  reflections  may  be  due  to  strong  impedance  contrasts,  such  as  those  in  a  salt 
dome  environment.  Single-station  deconvolutions  for  the  four  salt  shots  of  October  16,  1982 
recorded  at  the  two  SRO  stations  ANMO  and  KONO  are  also  obtained  with  similar  results. 

The  last  shot  in  the  sequence  of  October  16,  1982  appears  to  be  anomalous  with  respect  to  the 
three  earlier  shots,  as  evidenced  by  comparison  of  mean  amplitude  ratios  of  P  and  P-coda, 
spectral  ratio  slopes  of  P,  and  P/P-coda  spectral  slopes.  A  possible  reason  appears  to  be  altera¬ 
tion  of  shot  medium  due  to  close  proximity  of  an  earlier  shot. 

E.  ANALYSIS  OF  DATA  FROM  SOVIET  SALT  SHOTS  WITH  KNOWN  YIELDS  ,  / 

Teleseismic  P  arrivals  from  three  USSR  salt  explosions  with  known  yields  and  shot 
depths  and  recorded  at  the  LRSM  stations  NP-NT  and  RK-ON  are  analyzed.  Of  the  three 
Soviet  shots,  only  two  (22  April  1966  and  1  July  1968)  were  very  near  each  other  so  that  the 
propagation  path  effects  may  be  considered  to  be  the  same.  A  comparison  of  their  P  and  P- 
coda  spectral  ratios  at  both  stations  with  von  Seggem  and  Blandford  scaling  for  salt  indicates 
the  shot  medium  to  be  somewhat  "softer"  than  the  Salmon  salt.  The  deconvolution  method  of 
Shumway  and  Der  (1985)  is  applied  to  the  vertical  and  radial  component  records  at  NP-NT  of 
all  three  USSR  shots.  The  deconvolved  source  terms  for  only  the  two  deeper  (and  larger 
scaled  depth)  explosions  show  clear  evidence  of  P  and  pP  arrivals.  The  average  ratio  P/P-coda 
values  are  larger  for  explosions  with  larger  scaled  depths,  perhaps  because  an  overburied  shot 
is  likely  to  generate  less  near-source  coda  than  an  underburied  shot. 

security  ci ASSifiCAiiON  or  mis  page 
UNCLASSIFIED 


ANNUAL  REPORT 


STUDIES  IN  DECOUPLING 


TGAL-87-4 


MISCELLANEOUS  STUDIES  IN  DECOUPLING 
SUMMARY 

A.  COMPARISON  OF  SYNTHETIC  AND  OBSERVED  SEISMOGRAMS 

Synthetic  responses  for  plane  layered  model  applicable  for  paths  from  Salmon  to  dis¬ 
tances  of  16  and  32  km  (the  same  source-receiver  distances  as  for  stations  10S  and  20S, 
respectively)  are  computed  by  using  the  wavenumber  integration  method.  Several  arrivals  in 
the  synthetics  agree  well  with  the  observed  data  and  the  synthetic  seismograms  help  identify 
the  beginning  of  the  S  wavetrains  on  the  records  at  both  10S  and  20S.  Seismic  moment  of 
Salmon,  derived  by  a  comparison  of  the  peak  amplitudes  on  the  observed  and  theoretical 
seismograms,  agrees  remarkably  well  with  the  value  obtained  from  near-field  data. 

B.  ANALYSIS  OF  PROJECT  COWBOY  STRONG-MOTION  DATA 

Near-field  data  from  coupled  and  decoupled  explosions  in  salt,  acquired  under  Project 
Cowboy,  are  analyzed  by  obtaining  the  RVP  spectra  corrected  for  frequency-dependent 
geometrical  spreading.  Limited  data  suggest  strain  and  frequency-dependent  attenuation,  similar 
to  the  results  from  strong-motion  data  from  Salmon  and  Sterling. 

C.  EFFECT  OF  SCALED  DEPTH  ON  REGIONAL  PHASES 

Spectral  characteristics  of  the  regional  phases  Pn,  Pg,  and  Lg  are  studied  for  several  NTS 
explosions  covering  a  large  range  of  yields  and  scaled  depths.  A  new  parameter,  Effective 
Scaled  Depth  (ESD),  defined  as  the  ratio  of  shot  depth  to  the  estimated  elastic  radius,  is  intro¬ 
duced;  it  is  a  measure  of  scaled  depth  that  has  been  corrected  for  the  variation  in  overburden 
pressure  with  shot  depth.  Data  from  a  common  recording  station  indicate  the  mean  Pn/Lg  and 
Pg/Lg  spectral  slopes  to  be  most  dependent  on  ESD,  somewhat  less  on  scaled  depth,  and  the 
least  on  shot  depth.  The  spectral  shapes  of  Pn  and  Pg  are  significantly  more  dependent  on  ESD 
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and  scaled  depth  than  on  shot  depth.  The  enrichment  of  higher  frequencies  due  to  increased 
ESD  (or  scaled  depth)  is  the  largest  in  the  spectra  of  Pn,  somewhat  less  for  Pg,  and  the  least 
for  Lg.  The  low-frequency  asymptote  and  the  area  under  the  initial  pulse  of  deconvolved  Pn 
also  show  much  greater  dependence  on  scaled  depth  and  ESD  than  on  shot  depth. 

A  comparison  of  inter-event  spectral  ratios  of  Pn,  Pg,  and  Lg  for  shots  with  nearly  equal 
yields  but  different  ESD  or  scaled  depths  shows  significant  differences.  Similar  variations  are 
observed  if  PSRV  response  spectra  based  on  Pn  and  the  entire  seismogram  are  compared. 
These  results  imply  that  scaling  of  the  nuclear  seismic-source  function  may  vary  significantly 
from  one  regional  phase  to  another.  Therefore,  scaling  relationships  based  on  the  entire  seismo¬ 
gram  (e.g.  Mueller  and  Murphy,  1971),  generally  dominated  by  the  Lg  wavetrain,  may  not  be 
appropriate  for  Pn. 

D.  DECONVOLUTION  AND  SPECTRAL  RATIOS  FROM  SOVIET  SALT  SHOTS 

The  maximum-likelihood  multichannel  deconvolution  method  of  Shumway  and  Der 
(1985)  is  applied  to  teleseismic  P  waves  from  ten  salt  shots  in  the  Astrakhan  region  of  the 
USSR.  The  data  used  came  from  the  GBA  array  and  the  four  RSTN  stations,  RSON,  RSNT, 
RSNY,  and  RSSD.  The  deconvolved  source  functions  for  each  shot  show  clear  P  and  pP 
arrivals  separated  by  about  1  sec,  suggesting  large  shot  depths.  There  are  also  several  arrivals 
which  consistently  show  up  between  the  P  and  pP,  similar  to  those  observed  for  Salmon. 
These  earlier  reflections  may  be  due  to  strong  impedance  contrasts,  such  as  those  in  a  salt 
dome  environment.  Single-station  deconvolutions  for  the  four  salt  shots  of  October  16,  1932 
recorded  at  the  two  SRO  stations  ANMO  and  KONO  are  also  obtained  with  similar  results. 
The  last  shot  in  the  sequence  of  October  16,  1982  appears  to  be  anomalous  with  respect  to  the 
three  earlier  shots,  as  evidenced  by  comparison  of  mean  amplitude  ratios  of  P  and  P-coda, 
spectral  ratio  slopes  of  P,  and  P/P-coda  spectral  slopes.  A  possible  reason  appears  to  be 
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weakening  of  shot  medium  due  to  close  proximity  of  an  earlier  shot. 

E.  ANALYSIS  OF  DATA  FROM  SOVIET  SALT  SHOTS  WITH  KNOWN  YIELDS 

Teleseismic  P  arrivals  from  three  USSR  salt  explosions  with  known  yields  and  shot 
depths  and  recorded  at  the  LRSM  stations  NP-NT  and  RK-ON  are  analyzed.  Of  the  three 
Soviet  shots,  only  two  (22  April  1966  and  1  July  1968)  were  very  near  each  other  so  that  the 
propagation  path  effects  may  be  considered  to  be  the  same.  A  comparison  of  their  P  and  P- 
coda  spectral  ratios  at  both  stations  with  von  Seggem  and  Blandford  scaling  for  salt  indicates 
the  shot  medium  to  be  somewhat  "softer"  than  the  Salmon  salt.  The  deconvolution  method  of 
Shumway  and  Der  (1985)  is  applied  to  the  vertical  and  radial  component  records  at  NP-NT  of 
all  three  USSR  shots.  The  deconvolved  source  terms  for  only  the  two  deeper  (and  larger 
scaled  depth)  explosions  show  clear  evidence  of  P  and  pP  arrivals.  The  average  ratio  P/P-coda 
values  are  larger  for  explosions  with  larger  scaled  depths,  perhaps  because  an  overburied  shot 
is  likely  to  generate  less  near-source  coda  than  an  underburied  shot 
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LIST  OF  FIGURES 

Figure  No.  Title  Page 

A1  Synthetics  for  Salmon  at  16  km.  RAD,  TANG,  and  VERT  denote  radi-  6 

al,  tangential,  and  vertical  components,  respectively  whereas  VSS, 

VDS,  CLVD,  and  EXPL  represent  vertical  strike-slip,  vertical  dip-slip, 
compensated  vector-dipole,  and  explosion  sources,  respectively. 

A2  Synthetics  for  Sterling  at  16  km.  The  abbreviations  have  the  same  7 

meaning  as  in  Figure  Al. 

A3  Synthetics  for  Green’s  function  at  16  km.  The  abbreviations  have  the  8 

same  meaning  as  in  Figure  Al. 

A4  Synthetics  for  Salmon  at  32  km.  The  abbreviations  have  the  same  10 

meaning  as  in  Figure  Al. 

A5  Synthetics  for  Sterling  at  32  km.  The  abbreviations  have  the  same  11 

meaning  as  in  Figure  Al. 

A6  Synthetics  for  Green’s  function  at  32  km.  The  abbreviations  have  the  12 

same  meaning  as  in  Figure  Al. 

A7  Observed  and  synthetic  (explosion)  vertical  component  records  (with  13 

up  motion  representing  compression)  for  Salmon  at  the  same  time 
scale,  at  a  distance  of  16  km,  showing  good  agreement  between  ar¬ 
rivals  (such  as  those  marked  a,  b,  c,  and  d)  within  the  first  several 
seconds.  The  maximum  zero-to-peak  amplitude  values  are  given  at  the 
bottom  of  the  two  traces. 

A8  Observed  and  synthetic  (explosion)  vertical  component  records  (with  15 

up  motion  representing  compression)  for  Salmon  at  the  same  time 
scale,  at  a' distance  of  32  km,  showing  a  general  lack  of  agreement 
between  various  arrivals.  The  maximum  zero-to-peak  amplitude  values 
are  given  at  the  bottom  of  the  two  traces. 

B1  Record  Cow-30  (shown  at  bottom)  and  the  corresponding  RVP  spectra  18 

for  Shot  No.  9  (500  lbs,  coupled)  recorded  at  a  distance  of  15.4  m. 

The  lower  spectra  is  for  noise  preceding  the  signal  and  of  equal  win¬ 
dow  length. 

B2  Similar  to  and  for  the  same  shot  as  in  Figure  B1  for  Cow-29  and  a  dis-  19 

tance  of  24.3  m. 
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B3  Record  Cow-5  (shown  at  bottom)  and  the  corresponding  RVP  spectra  20 
for  Shot  No.  2  (20  lbs,  decoupled)  recorded  at  a  distance  of  6.3  m. 

The  lower  spectra  is  for  noise  preceding  the  signal  and  of  equal  win¬ 
dow  length. 


B4 


B5 


B6 

B7 


B8 

B9 
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Similar  to  and  for  the  same  shot  as  in  Figure  B3  for  Cow-3  and  a  dis-  21 
tance  of  1 1.0  m. 

Spectral  ratios,  corrected  for  noise,  based  on  RVP  spectra  derived  from  23 
Shot  9  at  distances  of  15.4  and  24.3  m  (upper  plot)  and  15.4  and  112.4 
m.  Points  for  which  S/N  power  ratio  is  less  than  2  are  not  plotted.  The 
dashed  line  shows  mean  (least-squares)  slope  over  the  frequency  range 
of  1  to  200  Hz. 

Similar  to  Figure  B5  for  Shot  2  (upper  plot)  and  Shot  8.  24 

Mean  values  of  100/Q  plotted  versus  average  source-receiver  distance  25 
for  the  frequency  range  of  1  to  100  Hz  for  both  coupled  and  decoupled 
shots.  The  vertical  lines  represent  error  bars  with  one  standard  devia¬ 
tion  whereas  the  horizontal  lines  denote  the  differential  source-receiver 
distance  over  which  mean  Q  has  been  computed. 

Similar  to  Figure  B7  but  for  the  frequency  range  of  1  to  200  Hz.  26 

Mean  attenuation  100/Q  over  the  frequency  range  of  1  to  25  Hz  versus  28 

source-receiver  distance,  R  from  Salmon  near-field  data  (after  Gupta  et 
al.,  1986).  The  horizontal  bars  denote  the  distance  range  and  the  verti¬ 
cal  lines  represent  error  bars  with  one  standard  deviation.  The  data  are 
consistent  with  a  R~1/2  decay  with  distance  (shown  by  continuous 
curve).  Note  that  results  in  Figures  B7  and  B8  show  variation  with  dis¬ 
tance  consistent  with  that  shown  above. 

Spectral  ratios  Duryea/Buteo  and  Cabriolet/Buteo,  corrected  for  noise,  33 
for  Pn,  Pg,  and  Lg,  recorded  at  KN-UT.  Points  for  which  S/N  power 
ratio  is  less  than  2  are  not  plotted.  The  mean  (least  squares)  slope  over 
the  frequency  range  1  to  5  Hz  (dashed  line)  is  indicated  (with  its  stan¬ 
dard  deviation  value  in  parentheses)  for  each  phase. 

Mean  Pn/Lg  spectral  slope  (per  Hz)  over  the  frequency  range  1  to  5  34 

Hz  versus  (a)  log  (shot  depth  in  meters),  (b)  scaled  depth,  and  (c)  log 
ESD  for  7  explosions,  numbered  as  in  Table  Cl.  The  vertical  lines 
represent  error  bars  with  one  standard  deviation.  The  correlation 
coefficient  and  mean  slope  value  of  the  least  squares  linear  trend 
(dashed  line)  is  indicated  for  each  plot. 
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C3  Similar  to  Figure  C2  for  the  spectral  ratio  Pg/Lg.  Note  that  the  mean  36 

slope  values  of  the  linear  trends  in  (a),  (b),  and  (c)  are  smaller  than  the 
corresponding  values  in  Figure  C2,  indicating  smaller  variation  in  the 
spectral  slope  Pg/Lg  as  compared  to  Pn/Lg. 

C4  Mean  spectral  slope  for  Pn  (corrected  for  noise,  yield,  and  t*)  over  the  38 

frequency  range  1  to  5  Hz  versus  (a)  log  (shot  depth  in  meters),  (b) 
scaled  depth,  and  (c)  log  ESD  for  7  explosions,  numbered  as  in  Table 
Cl.  The  vertical  lines  represent  error  bars  with  one  standard  deviation. 

The  correlation  coefficient  and  mean  slope  value  of  the  least  squares 
linear  trend  (dashed  line)  is  indicated  for  each  plot. 

C5  Similar  to  Figure  C4  for  Pg.  Note  that  the  mean  slope  values  of  the  39 

linear  trends  in  (a),  (b),  and  (c)  are  smaller  than  the  corresponding 
values  in  Figure  C4,  indicating  smaller  variation  in  the  spectra  of  Pg  as 
compared  to  Pn. 

C6  Log  (Aq/Y)  versus  (a)  shot  depth,  (b)  scaled  depth,  and  (c)  log  ESD  41 

for  7  explosions,  numbered  as  in  Table  Cl.  The  correlation  coefficient 
and  mean  slope  value  of  the  least  squares  linear  trend  (dashed  line)  is 
indicated  for  each  plot. 

C7  Similar  to  Figure  C6  for  log  (QJY).  Note  that  the  mean  slope  values  43 

of  the  linear  trends  in  (a),  (b),  and  (c)  are  smaller  than  the  correspond¬ 
ing  values  in  Figure  C6,  indicating  somewhat  smaller  variation  in 
log  (QJY)  as  compared  to  log  (AJY). 

C8  Plot  (log-log  scale)  of  the  area  under  the  initial  pulse  of  deconvolved  44 

Pn,  Q0,  versus  the  low-frequency  asymptote,  A0,  both  in  units  of  nm- 
sec.  The  near  equivalence  of  tl<e  two  measures  is  confirmed  by  the 
high  correlation  coefficient  value  and  slope  close  to  unity. 

C9  Comparison  of  Rex  (overburied),  scaled  Rex,  and  Schooner  (under-  47 

buried)  spectra  at  SE-6,  Las  Vegas  (Horizontal  Component),  after  Fig¬ 
ure  11,  Mueller  and  Murphy  (1971). 

CIO  PSRV  spectra  based  on  radial  component  data  at  KN-UT  for  (a)  Pn  48 

(6.4  sec  window)  and  (b)  total  (ALL)  seismogram  (102.4  sec  window) 
for  Rex  and  Schooner.  Only  spectral  points  with  S/N  power  ratio  of  at 
least  2  are  plotted. 

Cll  Similar  to  Figure  CIO  for  Buteo  (overburied)  and  Cabriolet  (under-  49 

buried). 
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Similar  to  Figure  CIO  for  Buteo  (overburied)  and  Cabriolet  (under-  51 
buried)  but  for  the  vertical  (Z)  components. 

Deconvolved  source  time  functions  for  the  four  Astrakhan  events  of  55 
October  16,  1982  recorded  at  GBA  (left)  and  RSTN  (center)  and 
"joint"  deconvolutions  (right)  based  on  the  deconvolved  waveforms 
from  both  GBA  and  RSTN.  The  arrows  denote  the  peaks  of  P  and  pP 
arrivals. 

Deconvolved  source  time  functions  for  the  six  Astrakhan  events  of  57 

September  24,  1983,  recorded  at  GBA.  The  arrows  denote  the  peaks 
of  P  and  pP  arrivals. 

Deconvolved  source  time  functions  for  the  ten  Astrakhan  events  of  Oc-  58 
tober  16,  1982,  and  September  24,  1983,  recorded  at  GBA. 

Deconvolved  source  time  functions  for  the  ten  Astrakhan  events  of  Oc-  59 
tober  16,  1982,  and  September  24,  1983,  recorded  at  the  three  RSTN 
stations  RSON,  RSNT,  and  RSNY. 

Deconvolved  source  time  functions  for  nine  Astrakhan  events  (four  on  60 
October  16,  1982,  and  five  on  September  24,  1983)  recorded  at  the 
four  RSTN  stations  RSON,  RSNT,  RSNY,  and  RSSD. 

Deconvolved  source  time  functions  for  the  four  Astrakhan  events  of  62 
October  16,  1982,  recorded  at  ANMO  and  KONO.  The  numbers 
denote  zero-to-peak  amplitudes  in  nanometers.  The  arrows  denote  the 
peaks  of  P  and  pP  arrivals. 

Inter-event  mean  spectral  amplitude  ratios  of  P,  in  log  units  and  over  63 
the  frequency  range  1  to  5  Hz,  for  4  RSTN  and  4  GBA  sensors.  The 
ratios  indicate  the  fourth  event  to  be,  on  the  average,  larger  than  the 
first  three  by  about  0.2  unit. 

Inter-event  mean  spectral  amplitude  ratios  of  P  coda,  in  log  units  and  65 
over  the  frequency  range  1  to  5  Hz,  for  4  RSTN  and  4  GBA  sensors. 

The  ratios  indicate  the  fourth  event  to  be,  on  the  average,  larger  than 
the  first  three  by  about  0.1  unit. 

Inter-event  mean  spectral  ratio  slopes  of  P,  in  log  units  per  Hz  and  66 
over  the  frequency  range  1  to  5  Hz,  for  4  RSTN  and  4  GBA  sensors. 

The  slope  values  indicate  the  fourth  event  to  have,  on  the  average, 
about  0.1  unit/Hz  lower  slope  than  the  first  three  shots. 
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DIO  Mean  P/P-coda  spectral  slopes,  in  log  units  per  Hz  and  over  the  ffe-  67 
quency  range  1  to  5  Hz,  for  the  four  shots  of  October  16,  1982  record¬ 
ed  at  the  four  RSTN  stations.  The  vertical  lines  represent  error  bars 
with  one  standard  deviation.  At  each  station,  the  fourth  shot  has  a 
lower  spectral  slope  value  than  for  the  first  three  shots. 

El  Vertical  component  records  of  four  salt  shots  at  the  LRSM  station  NP-  70 

NT.  A  denotes  the  epicentral  distance  in  degrees. 

E2  Records  of  the  same  four  salt  shots  as  in  Figure  El  at  the  LRSM  sta-  71 

don  RK-ON.  A  denotes  the  epicentral  distance  in  degrees. 

E3  Observed  spectral  ratios  for  P  and  P-coda  at  stations  NP-NT  and  RK-  72 

ON  and  theoretical  spectral  ratios 

E4  Deconvolved  source  functions  of  three  USSR  salt  shots  derived  from  74 

their  vertical  and  radial  component  records  at  NP-NT.  The  arrows 
denote  the  peaks  on  P  and  pP  arrivals.  The  letters  W,  h,  S.D.,  and 
P/P-coda  denote  yield  (kt),  shot  depth,  scaled  depth,  and  average  value 
. 1  over  the  frequency  range  of  0.5  to  3.0  Hz,  respectively. 
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INTRODUCTION 

One  of  the  most  difficult  problems  in  the  detection  and  discrimination  of  underground 
nuclear  explosions  is  that  posed  by  the  evasion  possibilities  of  decoupled  detonations.  Of  cru¬ 
cial  importance  is  the  spectral  shapes  of  the  decoupled  versus  coupled  or  tamped  explosions 
and  the  mechanism  responsible  for  the  differences  on  their  records  at  regional  and  teleseismic 
distances.  The  Russians  have  conducted  a  large  number  of  explosions  in  salt  under  their  Peace¬ 
ful  Nuclear  Explosions  (PNE)  program.  These  explosions  in  salt  are  still  continuing,  as  evi¬ 
denced  by  the  two  shots  of  April  19,  1987,  in  the  Urals  Mountains  region.  It  is  possible  that 
some  of  the  smaller  explosions  in  salt  were  decoupled.  It  is  therefore  important  to  understand 
the  seismic  waves  generated  by  Soviet  explosions  in  salt  With  these  objectives,  we  have  been 
carrying  out  a  number  of  miscellaneous  studies  related  to  decoupling.  Our  earlier  work  has 
already  been  reported  in  two  annual  reports  (Gupta  et  al.,  1986a,  and  Gupta  et  al.,  1986b). 

In  Gupta  et  al.  (1986a),  analyses  of  available  data  from  Salmon,  the  decoupled  shot  Ster¬ 
ling,  and  the  high-explosive  Sterling  Calibration  shot  recorded  at  USCGS  stations  10S  and 
20S  and  the  USGS  station  PL-MS  (at  shot-receiver  distances  of  16,  32,  and  27  km,  respec¬ 
tively)  were  carried  out  to  further  our  understanding  of  the  decoupling  of  explosions  in  salt 
The  spectra  and  spectral  ratios  were  obtained  for  both  P  and  the  surface-wave  group,  Lg,  for 
signal  windows  of  various  lengths.  The  results  confirmed  the  earlier  findings  of  Blandford  and 
Woolson  (1979)  indicating:  (1)  reduced  decoupling  at  higher  frequencies  for  both  P  and  Lg, 
and  (2)  relatively  greater  decoupling  for  Lg  than  for  P  at  high  frequencies.  At  each  station,  the 
spectral  ratio  P/Lg  was  found  to  have  a  more  positive  slope  for  the  decoupled  explosion  Ster¬ 
ling  than  for  the  tamped  shot  Salmon. 

An  attempt  was  made  to  determine  why  the  decoupling  ratio  is  different  for  P  and  the 
surface-wave  group,  Lg.  The  theoretical  work  of  Cisternas  (1964)  shows  that  for  Rayleigh 
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waves  produced  by  the  application  of  pressure  at  the  walls  of  spherical  cavities  whose  centers 
are  at  a  fixed  depth  from  the  free  surface,  the  larger  the  volume,  the  greater  the  enrichment  of 
high  frequencies.  In  other  words,  for  the  generation  of  surface  waves,  larger  volume  sources 
have  shallower  "effective"  source  depths.  For  a  fixed  source  function,  one  would  therefore 
expect  that,  in  comparison  to  Sterling,  a  considerably  larger  volume  around  the  shot  point  was 
involved  in  the  generation  of  surface  waves  by  Salmon  so  that  Lg  (Salmon)  should  be  richer  in 
high  frequencies  than  Lg  (Sterling),  as  observed.  Further  support  for  this  explanation  comes 
from  the  observed  increase  of  the  spectral  ratio  Pn/Lg  with  scaled  depth  for  7  Pahute  Mesa 
explosions  recorded  at  the  LRSM  station  KNUT. 

A  comparison  of  the  spectral  ratio  P/P-coda  at  teleseismic  distances  with  the  spectral  ratio 
Pn/Lg  at  regional  distances,  for  NTS  explosions,  shows  that  the  two  are  similar  in  their  depen¬ 
dence  on  "effective"  scaled  d_pth.  Variations  in  the  mean  spectral  slope  of  P/P-coda  may 
therefore  be  diagnostic  of  variations  in  the  spectral  ratio  Pn/Lg  or  P/Lg  similar  to  those 
observed  for  the  tamped  explosion  Salmon  and  the  decoupled  shot  Stirling.  Teleseismic  P 
arrivals  recorded  at  RKON  and  NPNT  from  4  Azgir  salt  explosions  covering  the  magnitude 
range  of  mb  =  4.7  to  5.5  were  analyzed.  The  spectral  ratios  P/P-coda  were  found  to  be  more 
positive  (indicating  P  to  be  richer  in  high-frequency  energy  as  compared  to  P-coda)  for  the  two 
smallest  shots  than  for  the  other  two  larger  explosions.  It  is  therefore  likely  that  the  two  smal¬ 
lest  shots  were  overburied  or  possibly  decoupled.  These  preliminary  results  suggest  a  method 
for  distinguishing  between  decoupled  and  normal  explosions  by  comparing  the  spectral  ratios 
Pn/Lg  on  regional  data  and/or  P/P-coda  on  teleseismic  data. 

A  summary  of  work  in  Gupta  et  al.  (1986b)  follows: 

An  understanding  of  the  frequency  dependence  of  decoupling  is  essential  for  monitoring 
nuclear  explosions  detonated  in  large  underground  cavities.  Analysis  of  available  data  from  the 
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nuclear  explosions  Salmon  (5.3  kt  in  salt)  and  Sterling  (0.38  kt  shot  in  the  cavity  created  by 
Salmon)  recorded  at  shot-receiver  distances  of  16,  32,  and  27  km,  confirmed  the  earlier 
findings  of  Blandford  and  Woolson  (1979)  indicating:  (1)  reduced  decoupling  of  Sterling  at 
higher  frequencies  for  both  P  and  "Lg,"  and  (2)  relatively  greater  decoupling  of  Sterling  for  Lg 
than  for  P  at  high  frequencies.  In  an  attempt  to  understand  why  the  frequency  dependence  of 
the  decoupling  ratio  is  different  for  P  and  the  Lg  group,  we  examined  the  spectral  content  of 
Pn  and  Lg  from  7  Pahute  Mesa,  NTS  explosions  covering  a  wide  range  of  scaled  depths  and 
recorded  at  a  common  station,  KN-UT.  The  spectral  ratio  Pn/Lg  varied  strongly  with  scaled 
depth,  mainly  due  to  significantly  greater  dependence  on  scaled  depth  of  the  spectra  of  Pn  than 
of  Lg.  These  results  are  in  agreement  with  those  from  Salmon  and  Sterling  with  scaled  depths 
of  475  and  1145  m/(kt)1/3,  respectively.  A  possible  reason  for  the  decoupling  at  higher  frequen¬ 
cies  to  be  different  for  P  and  Lg  is  therefore  the  greater  scaled  depth  of  decoupled  shots.  The 
larger  variability  of  Pn  spectra  than  those  of  Lg  can  be  due  to  an  explosion  (especially  those 
that  are  not  overburied)  approximating  a  pressurized  ellipsoid  cavity  with  major  axis  in  the 
vertical  direction  so  that  the  source  function  is  more  variable  for  down-going  waves  comprising 
Pn  than  for  others  making  up  Lg.  The  spectral  ratio  Pn/Lg  may  be  useful  in  discriminating 
between  decoupled  (or  overburied)  and  normal  shots. 

The  amplitude  and  spectra  of  teleseismic  and  regional  seismic  P  waves  depend  on  the 
material  properties  in  the  source  region.  The  effects  of  strain  and  frequency  dependent  attenua¬ 
tion  are  therefore  important  for  understanding  the  source  characteristics  of  coupled  versus 
decoupled  explosions.  Near-field  velocity  and  acceleration  data  from  the  Salmon  and  Sterling 
explosions  in  salt  were  analyzed  for  strain  and  frequency  dependent  attenuation.  Attenuation, 
parameterized  as  1/Q,  was  estimated  for  various  source-receiver  distance  ranges  and  frequency 
band  widths  by  using  two  methods:  (a)  mean  spectral  ratio  slope  and  (b)  average  pmplitude 
spectral  ratio.  Results  from  Salmon  indicate  that  (1)  in  the  frequency  range  of  1-25  Hz,  Q  is 
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abou*  5  to  10;  (2)  Q  increases  with  source-receiver  distance,  suggesting  lower  Q  for  larger 
strain  levels;  (3)  Q  in  the  frequency  range  of  25-50  Hz  is  substantially  higher  than  in  the  fre¬ 
quency  range  of  1-25  Hz,  or  Q  increases  with  frequency.  Analysis  of  data  from  Sterling  indi¬ 
cates  that  Q  is  frequency  independent  and  about  one  order  of  magnitude  larger  than  for  Sal¬ 
mon.  Therefore  the  near-source  attenuation  was  strain  and  frequency  dependent  and  hence 
non-linear  for  Salmon  out  to  a  range  of  600  m. 

Most  of  the  work  during  the  current  year  was  directed  towards  the  following  five  tasks: 

(A)  Comparison  of  synthetic  seismograms  for  Salmon  with  the  observed  data, 

(B)  Analysis  of  near-field  data  from  paired  cavity-detonated  and  coupled  (tamped)  explosions 
in  salt  recorded  under  Project  Cowboy, 

(C)  Effect  of  scaled  depth  on  regional  phases  by  analyzing  data  from  NTS  shots  with  known 
yields  and  depths, 

(D)  Deconvolution  and  spectral  ratios  of  P  waves  from  Soviet  salt  explosions  recorded  at 
arrays  and  stations  at  teleseismic  distances,  and 

(E)  Analysis  of  teleseismic  data  from  Russian  salt  explosions  of  known  yields. 
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A.  COMPARISON  OF  SYNTHETIC  AND  OBSERVED  SEISMOGRAMS 

Theoretical  Seismograms 

Using  wavenumber  integration,  responses  for  plane  layered  earth  model  appropriate  for 
paths  within  32  km  from  Salmon  and  Sterling  were  computed  (C.  A.  Langston,  written  com¬ 
munication).  The  sampling  rate  is  20/sec  so  that  the  Nyquist  frequency  is  10  Hz.  The  earth 
model  (Table  Al)  was  based  on  Springer’s  (1966)  Model  D  which  is  essentially  derived  from 
Warren  et  al.’s  (1966)  refraction  results  in  the  area.  The  attenuation  parameter  Q  for  P  waves, 
QP,  is  assumed  to  be  twice  that  for  S  waves,  i.e.  QP  =  2  Qs. 


TABLE  Al 


EARTH  MODEL  USED  FOR  SYNTHETIC  SEISMOGRAMS 


.  P-wave  S-wave  ~ 

Thlckness  Velocity  QP  Velocity  Qs  DenS,,y 

(km)  (km/sec)  (km/sec)  (gm/cc) 


The  displacement  field  from  the  simplest  of  sources  -  namely,  the  unidirectional  unit 
impulse,  which  is  localized  precisely  in  both  space  and  time,  is  the  elastodynamic  Green  func¬ 
tion  (Aki  and  Richards,  1980).  Displacement  from  realistic  source  models  can  therefore  be 
obtained  by  synthesizing  appropriate  Green’s  functions.  The  synthetics  for  source-receiver  dis¬ 
tance  of  16  km  (the  same  as  for  records  at  station  10S)  are  shown  in  Figures  Al,  A2,  and  A3 
for  Salmon,  Sterling,  and  Green’s  function.  A  seismic  moment  of  1  x  1020  dyne-cm  is  assumed 


for  each  source.  Salmon  and  Sterling  have  Brune  source  functions  with  corner  frequencies  of 
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Figure  A3.  Synthetics  for  Green’s  function  at  16  km.  The  abbreviations  have  the  same 
meaning  as  in  Figure  Al. 
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j 

i 

j  4  and  10  Hz,  respectively.  The  synthetics  represent  ground  velocity  and  the  peak  amplitudes 

I 

are  normalized  with  the  scale  factor  (cm/sec)  shown  to  the  right.  The  top  four  traces  are  the 
radial  (RAD)  components  for  sources  needed  to  construct  an  arbitrary  moment  tensor  source; 

j 

|  VSS,  VDS,  CLVD,  and  EXPL  denote  vertical  strike-slip,  vertical  dip-slip,  compensated  linear 

I 

vector  dipole,  and  isotropic  sources,  respectively.  The  middle  two  traces  in  each  figure 
represent  the  tangential  (TANG)  components  whereas  the  last  four  show  the  vertical  (VERT) 
component  ground  velocity.  Conventions  for  the  coordinate  system  are  positive  radial  away 
from  the  source,  vertical  positive  downv/ard,  and  tangential  positive  clockwise  about  the  source 
(looking  from  above).  Results  for  the  same  four  sources  at  distance  of  32  km  (the  same  as  for 

t 

| 

|  the  20S  records)  are  shown  in  Figures  A4,  A5,  and  A6. 

i 

!  Comparison  with  Observed  Records  at  Stations  10S  and  20S 

f 

The  instrument  response  for  the  recording  at  10S  is  nearly  flat  in  ground  velocity  over  a 
wide  range  of  frequencies  (see  Springer  et  al.,  1968,  Figure  3).  A  comparison  of  the  observed 

! 

vertical  component  data  from  Salmon  at  one  of  the  detectors  at  station  10S  and  the  synthetic 
!  for  VERT  EXPL  (Figure  Al),  after  applying  the  instrumental  response,  is  shown  in  Figure  A7. 

Considering  the  simplicity  of  the  assumed  theoretical  model,  agreement  between  the  theoretical 

t 

and  observed  records  is  good  for  the  first  several  seconds  after  the  first  P  arrival.  Synthetics  for 

I 

the  TANG  VDS  source  (Figure  Al),  which  does  not  generate  direct  P  but  shows  a  distinct  first 

S  arrival,  indicate  an  S-P  time  of  about  3.5  sec.  The  largest  arrival  on  the  observed  record  is 

i 

the  low-frequency  arrival  starting  about  1.7  sec  after  the  onset  of  P.  This  maxima  in  the  wave 
trains  was  considered  to  be  the  S-wave  group  by  Springer  et  al.  (1968).  The  synthetics,  on  the 
other  hand,  suggest  that  this  arrival,  being  too  early  to  be  S  wave,  is  due  to  P  waves  rever- 

( 

berating  within  the  low  velocity  sediments.  Furthermore,  this  arrival  is  too  early  to  be  the 

direct  S  on  the  basis  of  the  available  three-component  data  from  Salmon  and  Sterling  at 

I 
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comparable  source-receiver  distances  in  Borcherdt  et  al.  (1967).  A  distinct  arrival,  about  3.3 
sec  after  the  first  P,  is  therefore  identified  as  the  beginning  of  the  S  wavetrain  (denoted  by  S  in 
Figure  A7a).  Beyond  the  arrival  of  S  waves,  the  synthetic  record  shows  lot  more  motion  than 
the  observed  record.  This  indicates  that  there  may  be  too  many  interfaces  within  the  model  of 
sediments  and/or  that  the  assumed  Q  is  too  high  so  later  arrivals  including  the  surface  waves 
are  not  attenuated  enough  at  high  frequencies. 

The  observed  and  corresponding  synthetic  vertical  component  records  from  station  20S 
are  shown  in  Figure  A8.  Agreement  between  the  observed  and  theoretical  seismograms  is  not 
as  good  as  for  the  data  from  station  10S.  Synthetics  (Figure  A4)  indicate  an  S-P  time  of  about 
5.8  sec  which  agrees  well  with  the  distinct  arrival  (denoted  by  S  in  Figure  A8a)  separated  from 
the  first  P  by  nearly  5.5  sec. 

Both  Figures  A7  and  A8  indicate  substantially  larger  S  waves  in  the  observed  data  than 
in  the  synthetics  for  an  explosion  source.  This  means  that  the  source  function  for  Salmon  was 
considerably  more  complicated  than  an  isotropic  point  source.  A  comparison  of  the  zero-to- 
peak  amplitude  values  in  Figures  A7  and  A8  suggest  the  ratio  observed/synthetic  of  about  104 
and  163,  respectively.  Since  the  synthetics  are  based  on  a  seismic  moment  of  1  x  1020  dyne- 
cm,  the  observed  data  from  10S  and  20S  indicate  seismic  moments  of  about  1.0  x  1022  and 
1.6  x  1022  dyne-cm,  respectively.  These  values  agree  remarkably  well  with  the  Salmon  seismic 
moment  value  of  2  x  1022  dyne-cm  derived  from  close-in  data  (Aki  et  al.,  1974). 

Conclusion 

Synthetic  seismograms  for  Salmon  show  considerably  less  S  waves  than  the  observed 
data  at  both  10S  and  20S.  This  means  that  the  source  function  was  significantly  more  compli¬ 
cated  than  an  isotropic  point  explosion  source.  The  synthetics  were  helpful  in  the  correct 
identification  of  the  beginning  of  the  S  wavetrains. 
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Figure  A8.  Observed  and  synthetic  (explosion)  vertical  component  records  (with  up  motion 
representing  compression)  for  Salmon  at  the  same  time  scale,  at  a  distance  of  32 
km  showing  a  general  lack  of  agreement  between  various  arrivals.  The  maximum 
zero-to-peak  amplitude  values  are  given  at  the  bottom  of  the  two  traces. 
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B.  ANALYSIS  OF  PROJECT  COWBOY  STRONG-MOTION  DATA 

Introduction 

Project  Cowboy  was  an  experiment  designed  to  determine  to  what  extent  underground 
explosions  could  be  effectively  concealed  simply  by  firing  the  explosion  in  a  large  cavity. 
Comparative  measurements  of  ground  motion  were  obtained  from  tamped  charges  of  high 
explosives  and  high-explosive  charges  of  similar  yield  fired  in  large  cavities  in  halite.  The 
tamped  explosions  are  referred  to  as  coupled  shots  and  the  explosions  in  cavities  as  decoupled 
shots.  The  experiment  was  conducted  at  depths  near  800  ft  in  a  salt  mine  near  Winnfield, 
Louisiana.  Details  of  the  experiment  have  been  provided  by  Murphey  (1961),  Adams  and  Allen 
(1960),  and  many  others.  Analysis  of  strong-motion  data  recorded  at  distances  less  than  600  ft 
was  first  performed  by  Murphey  (1961).  His  study  was,  however,  based  on  peak  measurements 
and  did  not  include  spectral  analyses. 

The  Project  Cowboy  data  were  digitized  at  5000  samples/sec  and  made  available  to  us 
(Neil  K.  Perl,  written  communication).  The  measuring  devices  were  two  different  types  of 
velocity  gages  with  different  instrumental  responses  (Thombrough  et  al.,  1960).  An  attempt 
was  made  to  analyze  as  much  data  as  possible  but  a  large  number  of  recordings  suffered  from 
problems  such  as  clipping  and/or  poor  S/N  ratio.  Many  records  of  decoupled  shots  showed 
high  frequency  signals,  probably  due  to  ringing  of  the  canister  containing  the  gages  (Murphey, 
1961).  Useful  data  could  therefore  be  obtained  from  only  a  rather  small  number  (13)  of  record¬ 
ings  as  listed  in  Table  Bl. 
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TABLE  B1 

PROJECT  COWBOY  DATA  USED  IN  THIS  STUDY 


Record  No. 

Shot  No. 

Yield  (lb) 

Coupled  (C) 

or  Distance  (m) 

Decoupled  (D) 

Cow-3 

2 

20 

D 

11.0 

Cow-5 

2 

20 

D 

6.3 

Cow- 12 

6 

200 

D 

24.3 

Cow-14 

6 

200 

D 

61.0 

Cow-15 

6 

200 

D 

24.4 

Cow-23 

8 

477 

D 

24.3 

Cow-26 

8 

477 

D 

61.0 

Cow-29 

9 

500 

C 

24.3 

Cow-30 

9 

500 

C 

15.4 

Cow-32 

9 

500 

C 

112.4 

Cow-72 

17 

200 

C 

45.7 

Cow-73 

17 

200 

C 

50.9 

Cow-76 

17 

200 

C 

83.5 

Attenuation  from  Coupled  and  Decoupled  Shots 

By  analyzing  near-field  velocity  and  acceleration  data  from  the  Salmon  and  Sterling 
explosions  in  salt,  Gupta  et  al.  (1986b)  obtained  evidence  for  strain  and  frequency  dependent 
attenuation.  An  attempt  was  made  to  carry  out  similar  analysis  by  obtaining  spectral  ratios 
between  various  receiver  locations  for  both  coupled  (tamped)  and  decoupled  shots.  Before  the 
spectral  ratios  were  taken,  it  was  of  course  necessary  to  correct  the  observed  records  from 
identical  instruments  for  the  distance  and  frequency-dependent  geometrical  spreading.  An  algo¬ 
rithm  for  applying  this  correction  (e.g.  Larson,  1982)  and  thereby  deriving  spectra  that 
correspond  to  spectra  of  the  Reduced  Velocity  Potential  (RVP)  was  developed.  Four  examples 
are  shown  in  Figures  Bl,  B2,  B3,  and  B4  in  which  the  first  two  are  from  a  coupled  shot  (Shot 
No.  9,  500  lbs)  recorded  at  distances  of  15.4  and  24.3  m,  respectively.  The  other  two  figures 


Teledyne  Geotech 


17 


November  1987 


AMPLITUDE  CCM**3/SEC]  *  SQP 


0.00  40.  QQ  80.00  120.00 


Wed  Aug  12  11  :0'5  =4 1 
COWBOY244 
SEIS  #  244S 
WINDOW  =  1024 
SIGNAL  1072  1024 
NOISE  148  1024 

30. oox  cos  r  A  PE  II 

TAPER  RIGHT  END  ONLY 

DTRND2  1024  POINTS 

NO  INST  CORRECTION 

NO  SMOOTHING 

RVP»  4400.  15.4 

NO  DECIMATION 

+  +  c  ou;3Q 

+  + 

+  2.3-2 -V 

+  +  +  +  + 

++  COWBOY  9  SHOT 


160.00  200.00 


FREQUENCY  Hz 


O0O8  SEC  1 2 Mo  92  CM/SEC 

Figure  Bl.  Record  Cow-30  (shown  at  bottom)  and  the  corresponding  RVP  spectra  for  Shot 
No.  9  (500  lbs,  coupled)  recorded  at  a  distance  of  15.4  m.  The  lower  spectra  is 
for  noise  preceding  the  signal  and  of  equal  window  length. 
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Figure  B2.  Similar  to  and  for  the  same  shot  as  in  Figure  B1  for  Cow-29  and  a  distance  of 
24.3  m. 
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Figure  B3.  Record  Cow-5  (shown  at  bottom)  and  the  corresponding  RVP  spectra  for  Shot 
No.  2  (20  lbs,  decoupled)  recorded  at  a  distance  of  6.3  m.  The  lower  spectra  is 
for  noise  preceding  the  signal  and  of  equal  window  length. 
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Figure  B4.  Similar  to  and  for  the  same  shot  as  in  Figure  B3  for  Cow-3  and  a  distance  of 
1 1.0  m. 
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are  from  a  decoupled  shot  (Shot  No.  2,  20  lbs  in  a  cavity  of  12  ft  diameter)  recorded  at  dis¬ 
tances  of  6.3  and  11.0  m,  respectively.  The  raw  trace,  with  its  peak  amplitude  value,  is  shown 
at  the  bottom  of  each  spectrum.  Spectral  ratios  from  the  two  spectra  in  Figures  B1  and  B2, 
corrected  for  noise  and  considering  only  those  points  with  S/N  power  ratio  of  at  least  2,  are 
shown  in  Figure  B5  (upper  plot),  in  which  the  dashed  line  shows  mean  (least  squares)  slope 
over  the  frequency  range  of  1  to  200  Hz.  Spectral  ratio  from  the  spectrum  of  Figure  B1  and 
another  at  a  distance  of  112.4  m  are  also  shown  in  Figure  B5  (bottom  plot).  Similar  spectral 
ratios  from  two  decoupled  shots  are  shown  in  Figure  B6,  in  which  the  dashed  lines  again 
represent  the  mean  slopes  over  the  same  frequency  range. 

The  mean  spectral  slopes  can  be  used  to  estimate  the  mean  Q,  Q,  over  a  given  frequency 
range.  The  mean  t*,  t*,  is  given  by  (e.g.  Der  and  Lees,  1985) 

t*  =  -  — .  d('nA)  (Bl) 

n  df 

where  A  and  f  represent  amplitude  and  frequency,  respectively.  If  the  amplitude  ratio  (in  log10 
units)  from  two  receivers  at  source-receiver  distances  Rj  and  R2  is  plotted  versus  f  (as  in  Fig¬ 
ures  B5  and  B6)  and  the  mean  slope  per  Hz  is  s,  then 

t*  =  0.733  s  (B2) 

and 

Q  =  -4r  (B3) 

u  t* 

where  A  =  R2  -  Rj  and  u  is  the  wave  velocity,  about  4.4  km/sec  for  salt  (Murphey,  1961).  In 
this  manner,  data  from  decoupled  shots  (Table  Bl)  were  used  to  obtain  the  four  spectral  ratios 
Cow-5/Cow-3,  Cow- 12/Cow- 14,  Cow- 15/Cow- 14,  and  Cow-23/Cow-26,  and  the  corresponding 
Q  values  were  estimated  for  frequency  ranges  of  1  to  100  Hz  and  1  to  200  Hz.  Similarly,  the 
spectral  ratios  Cow-30/Cow-29,  Cow-30/Cow-29,  Cow-72/Cow-76,  and  Cow-73/Cow-76  were 
used  to  derive  Q  values  from  coupled  shots.  The  results  are  shown  in  Figures  B7  and  B8,  in 
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Figure  B6.  Similar  to  Figure  B5  for  Shot  2 
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Figure  B8.  Similar  to  Figure  B7  but  for  the  frequency  range  of  1  to  200  Hz. 
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which  the  mean  values  of  100/Q  are  plotted  versus  source-receiver  distance;  the  vertical  lines 
represent  error  bars  with  one  standard  deviation,  whereas  the  horizontal  lines  denote  A.  The 
100/Q  values  are  plotted  at  the  mid-point  of  A  for  both  coupled  and  decoupled  shots. 

In  Gupta  et  al.’s  (1986b)  analysis  of  strong-motion  data  from  Salmon  and  Sterling  to 
obtain  evidence  for  strain  and  frequency-dependent  attenuation,  average  1/Q  was  measured  for 
each  distance  range  bracketed  by  the  sensor  distances.  Figure  B9,  reproduced  from  their  study, 
shows  the  mean  100/Q  (a  measure  of  attenuation)  from  Salmon  to  vary  with  strain,  approxi¬ 
mately  as  R-1/2.  The  data  in  Figures  B7  and  B8  are  too  limited  to  draw  any  firm  conclusions 
but  appear  to  be  similar  to  those  in  Figure  B9,  with  attenuation  decreasing  with  R  for  both 
coupled  and  decoupled  shots.  At  comparable  distances,  the  mean  100/Q  values  for  coupled 
shots  seem  to  be  larger  than  those  for  decoupled  shots,  probably  because  of  larger  strain  asso¬ 
ciated  with  coupled  shots.  Furthermore,  the  mean  attenuation  values  for  the  frequency  range  of 
1  to  100  Hz  are  larger  than  those  for  the  frequency  range  of  1  to  200  Hz,  perhaps  due  to  Q 
increasing  with  frequency,  as  also  observed  in  the  Salmon  data  (Gupta  et  al.,  1986b). 

Conclusion 

Data  from  both  coupled  and  decoupled  shots  appear  to  support  our  earlier  results  (Gupta 
et  al.,  1986b)  suggesting  strain  and  frequency-dependent  attenuation.  Note,  however,  that  the 
results  here  are  not  only  based  on  very  limited  data  but  also  suffer  from  inaccuracies  due  to 
spectral  slopes  being  measured  over  rather  small  distance  range  (from  about  5  to  100m). 
Another  source  of  error  is  the  presence  of  prominent  secondary  arrivals  (such  as  reflection 
from  the  mine  roof)  within  the  time  window  used  for  spectral  analysis,  especially  at  larger 
source-receiver  distances. 
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Figure  B9.  Mean  attenuation  100/Q  over  the  frequency  range  of  1  to  25  Hz  versus 
source-receiver  distance,  R  from  Salmon  near-field  data  (after  Gupta  et  al.,  1986b). 
The  horizontal  bars  denote  the  distance  range  and  the  vertical  lines  represent  error 
bars  with  one  standard  deviation.  The  data  are  consistent  with  a  R~,/2  decay  with 
distance  (shown  by  continuous  curve).  Note  that  results  in  Figures  B7  and  B8 
show  variation  with  distance  consistent  with  that  shown  above. 
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C.  EFFECT  OF  SCALED  DEPTH  ON  REGIONAL  PHASES 

Introduction 

Numerous  observational  studies  of  regional  phases  (e.g.  Blandford,  1981,  1982;  Barker  et 
al.,  1981;  Pomeroy  et  al.,  1982;  Gupta  et  al.,  1982;  Gupta  and  Blandford,  1983;  Der  et  al., 
1984;  Glaser  et  al.,  1986;  Gupta  and  McLaughlin,  1987)  indicate  large  variations  in  amplitude 
and  spectra  due  to  factors  such  as  site  effects,  sensitivity  to  near-source  environment,  and  pro¬ 
pagation  path  effects.  Extreme  variability  of  regional  phases  from  NTS  explosions  at  essen¬ 
tially  the  same  location  has  been  noted  by  Springer  and  Denny  (1976,  especially  their  Figure  4) 
and  Gupta  et  al.  (1984,  see  their  Figure  11).  Theoretical  studies  based  on  synthetic  seismo- 
grams  of  regional  phases  (e.g.  Bouchon,  1982;  Banda  et  al.,  1982;  Campillo  et  al.,  1984; 
Rivers,  1984,  Bennett  et  al.,  1987)  also  suggest  that  the  excitation  and  propagation  of  regional 
phases  can  be  highly  variable.  In  this  study,  we  have  tried  to  isolate  the  effects  of  source  and 
near-source  parameters,  especially  depth  and  yield,  by  comparing  regional  phases  from  shots 
within  the  Pahute  Mesa  region  and  recorded  at  a  common  station. 

Effect  of  Scaled  Depth  on  Spectra  of  Regional  Phases 

Scaled  depth  for  nuclear  explosions  is  generally  defined  as  depth(m)/YI/3  where  Y  is 
yield  in  kton.  Several  studies  have  indicated  that  scaled  depth  is  an  important  parameter 

i 

directly  influencing  the  spectra  and  amplitudes  of  the  explosion-generated  body  and  surface 
waves  (e.g.  Mueller  and  Murphy,  1971;  Rulev,  1965;  and  Kisslinger,  1963  for  a  review  of 
early  work).  The  elastic  radius,  r,  varies  as  Y1/3  so  that  scaled  depth  varies  as  h/r  where  h  is 
the  shot  depth.  The  elastic  radius,  r,  should  also  vary  inversely  as  rock  strength  since  the  non¬ 
linear  zone  extends  to  the  point  where  the  explosive  hoop  stress  equals  the  tensile  strength  of 

| 

the  surrounding  rock.  This  results  in  larger  elastic  radii  for  shots  in  lower-strength  media 
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(Evemden  et  al.,  1986).  Thus  r  would  depend  not  only  on  Y  but  also  on  medium  properties 
such  as  density  and  overburden  pressure  which  are  related  to  the  rock  strength.  For  example, 
Mueller  and  Murphy  (1971)  empirically  determined  that  for  tuff-rhyolite  detonations, 


r  _  2000 

Y1/3  (ph)I/24 


(Cl) 


where  p  is  the  mean  overburden  (average  shot  point  to  surface)  density  in  gm/cc  and  h  is  in 
meters.  We  shall  refer  to  the  quantity  h/r  as  "Effective  Scaled  Depth"  or  ESD  so  that,  for 
tuff-rhyolite  shots,  ESD  may  be  approximated  by 


ESD  = 


(C2) 


The  first  term  on  the  right  hand  side  of  the  above  equation  is  the  scaled  depth,  whereas  the 
second  term  may  be  considered  a  "correction"  for  the  variation  in  overburden  pressure  as  the 
shot  depth  varies. 

If  explosions  at  considerably  different  depths  and  in  significantly  different  media  are  to 
be  compared,  physical  considerations  would  suggest  the  spectral  content  of  the  radiated  seismic 
waves  to  be  controlled  by  the  dimensionless  quantity  h/r  or  ESD,  rather  than  by  the  scaled 
depth  h/Y1/3.  Note  that  Cistemas’  (1964)  study  of  the  radiation  of  elastic  waves  from  a  spheri¬ 
cal  cavity  in  a  homogeneous  half  space  also  showed  the  surface  displacement  to  be  dependent 
on  the  ratio  h/r. 


In  order  to  study  the  effect  of  variations  in  scaled  depth,  we  shall  now  describe  some 
observational  results  comparing  the  spectral  contents  of  Pn,  Pg,  and  Lg  for  explosions  with 
large  differences  in  their  mot  depth,  scaled  depth  and  ESD  values.  The  Pahute  Mesa  explo¬ 
sions  Buteo,  Rex,  Scotch,  and  Benham,  covering  the  yield  range  of  0.7  to  1100  kt,  provide 
unique  data  for  a  study  of  amplitude-yield  relations  at  regional  distances  (Blandford,  1976). 
Regional  phases  from  these  four  and  three  additional  explosions  listed  in  Table  Cl  were  well 
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Name 


Date 


Yield+ 


ESD 


Shot 

Medium 


Shot 

Depth 


Depth  of 
Water  Table 


Scaled 

Depth# 


(m) 


(m) 


(kt) 


Buteo 

Rex 

Duryea 

Scotch 

Cabriolet* 

Schooner* 

Benham 


12  May  65 
24  Feb  66 
14  Apr  66 
23  May  67 
26  Jan  68 
8  Dec  68 
19  Dec  68 


tuff 

tuff 

rhyolite 

tuff 

rhyolite 

tuff 

tuff 


696 

671 

544 

977 

52 

111 

1402 


660 

642 

662 

672 

488 

274 

641 


0.7 

784 

8.049 

16 

266 

2.682 

65 

135 

1.222 

150 

184 

2.262 

2.3 

39.4 

0.123 

35 

33.9 

0.146 

1100 

136 

2.053 
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To  isolate  source  effects,  it  is  important  to  use  common  source-receiver  paths  because  of 
the  extreme  sensitivity  to  propagation  paths  of  the  regional  phases,  especially  Pn  (Langston, 
1982)  and  Lg  (Gupta  and  Blandford,  1983).  The  seven  explosions  in  Table  Cl  are  all  from 
Pahute  Mesa  so  that  the  source-receiver  paths  to  KN-UT  should  not  be  much  different,  and  the 
results  should  be  nearly  free  from  propagation  path  effects.  The  window  lengths  used  for 
spectral  analysis  of  Pn,  Pg,  and  Lg  on  the  vertical  component  records  were  6.4,  12.8,  and  12.8 
sec,  respectively,  and  a  10%  cosine  taper  was  used  along  with  3,  5,  and  5  point  smoothing  of 
spectra  for  Pn,  Pg,  and  Lg  windows,  respectively.  Out  of  these  7  shots,  Buteo  and  Duryea 
(with  their  yields  differing  by  a  factor  of  about  100)  were  detonates  in  the  same  hole  so  that 
the  source-receiver  paths  are  almost  identical.  A  comparison  of  the  relative  variations  in  spectra 
of  the  three  regional  phases  for  various  explosions  can  be  made  by  an  examination  of  inter¬ 
event  spectral  ratios.  Two  such  ratios  Duryea/Buteo  and  Cabriolet/Buteo  for  Pn,  Pg,  and  Lg, 
each  corrected  for  noise,  are  shown  in  Figure  Cl.  Cabriolet  and  Buteo  differ  in  their  yield 
values  by  a  factor  of  only  about  3.  For  both  sets  of  inter-event  spectral  ratios,  the  mean  spec¬ 
tral  slopes  over  the  frequency  range  of  1  to  5  Hz  are  the  largest  for  Pn,  the  least  for  Lg,  and 
intermediate  between  the  two  for  Pg.  It  seems  therefore  that  the  spectra  of  Pn  are  the  most 
variable  and  those  of  Lg  the  least  variable. 

We  next  obtained  the  spectral  ratios  Pn/Lg  at  the  common  station  KN-UT  for  the  seven 
explosions  listed  in  Table  Cl.  Results  are  shown  in  Figure  C2,  in  which  the  mean  spectral 
slopes  of  Pn/Lg  over  the  frequency  range  of  1  to  5  Hz  are  plotted  versus  (a)  shot  depth,  (b) 
scaled  depth,  and  (c)  ESD,  all  in  log  units.  The  three  plots  show  linear  trends  with  correlation 
coefficient  values  of  0.726,  0.811,  and  0.835,  respectively.  Assuming  the  validity  of  a  linear 
relationship,  it  seems  that  the  spectral  ratio  Pn/Lg  is  influenced  the  most  by  ESD,  somewhat 
less  by  scaled  depth,  and  the  least  by  shot  depth.  The  values  of  mean  slope  of  the  linear  trends 
in  the  three  plots  in  Figure  C2  are  also  indicated. 
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SLOPE  =  -0.134  (0.020) 


Figure  Cl.  Spectral  ratios  Duryea/Buteo  and  Cabriolet/Buteo,  corrected  for  noise,  for  Pn, 
Pg,  and  Lg,  recorded  at  KN-UT.  Points  for  which  S/N  power  ratio  is  less  than  2 
are  not  plotted.  The  mean  (least  squares)  slope  over  the  frequency  range  1  to  5 
Hz  (dashed  line)  is  indicated  (with  its  standard  deviation  value  in  parentheses)  for 
each  phase. 
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The  spectral  ratios  Pg/Lg  were  obtained  for  6  shots;  Pg  for  Rex  was  clipped  and  could 
not  be  used.  Plots  based  on  Pg/Lg  but  similar  to  those  in  Figure  C2  are  shown  in  Figure  C3. 
The  three  linear  trends  with  correlation  coefficients  of  0.849,  0.853,  and  0.909,  respectively, 
imply  that  the  spectral  ratio  Pg/Lg  is  also  influenced  the  most  by  ESD,  somewhat  less  by 
scaled  depth  and  the  least  by  shot  depth.  The  mean  slope  values  of  the  linear  trends  in  the 
three  plots,  indicated  in  Figure  C3,  are  all  considerably  smaller  than  the  corresponding  values 
for  Pn/Lg  in  Figure  C2.  This  means  that,  compared  to  Pn/Lg,  Pg/Lg  has  a  somewhat  weaker 
dependence  on  shot  depth,  scaled  depth,  and  ESD. 

We  also  compared  the  spectral  content  of  individual  Pn  and  Pg  windows  after  correcting 
for  the  large  differences  in  the  explosion  yields.  For  the  latter,  we  made  use  of  von  Seggem 
and  Blandford’s  (1972)  source  scaling,  in  which  the  source  function  for  a  given  yield  Y  is 
expressed  as 

T(°o)  (C4) 

[«B/k)2  +  1  J 

where  (0  =  2  tc  f,  'F(°°)  is  a  coupling  term  independent  of  frequency, 

k  =  ko  (5/Y),/3  (C5) 

and 

A  =  2B  +  1  (C6) 

where  B  and  kQ  are  medium-dependent  constants.  In  this  model,  shot  depth  is  not  an  indepen¬ 
dent  parameter,  but  the  source  function  can  be  varied  by  assigning  different  values  to  the  rise 
time  and  overshoot  parameters  ko  and  B,  respectively.  Considering  first  the  Pn  windows  for  the 
various  explosions,  the  spectra  were  corrected  for  the  instrument  response  and  anelastic 
attenuation  by  assuming  t*  =  0.1  sec,  the  value  used  by  Blandford  (1976).  A  correction  for  the 
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differences  in  yield  was  applied  by  using  the  source  function  in  equation  (C4)  to  correspond 
to  a  tuff  reduced  displacement  potential  (RDP),  i.e.  ko  =  12.0  and  B  =  0.  With  these  three 
corrections,  the  Pn  spectra  become  flat  if  the  shot  depth  and  medium  did  not  influence  the 
spectral  shapes.  Most  of  the  corrected  spectra  were  not  flat  and  their  mean  spectral  slopes  indi¬ 
cated  nearly  linear  trends  (the  results  for  Duryea  and  Buteo  can  be  seen  in  Gupta  et  al.,  1986b, 

Figure  18).  This  means  that  ordinary  cube-root  scaling  (i.e.  with  fixed  values  of  parameters  ko 
and  B)  fails  to  model  the  observed  Pn  spectra  of  explosions  with  large  variations  in  yield  and 
shot  depth.  Plots  of  the  derived  mean  spectral  slopes  (over  the  frequency  range  of  1  to  5  Hz) 
versus  the  shot  depth,  scaled  depth,  and  ESD  (all  in  log  units)  are  shown  in  Figure  C4.  The 
least  squares  linear  trends  (dashed  lines)  have  correlation  coefficients  of  0.827,  0.970,  and 
0.979,  respectively,  suggesting  that  the  Pn  spectra  are  influenced  the  most  by  ESD,  somewhat 
less  by  scaled  depth,  and  the  least  by  shot  depth.  ] 

I 

The  same  corrections  (i.e.  t*  =  0.1  sec  and  tuff  RDP)  were  also  applied  to  the  Pg  spectra 
and  the  mean  spectral  slopes  over  the  same  frequency  range  computed.  Plots  of  these  slopes 
versus  the  shot  depth,  scaled  depth,  and  ESD  are  shown  in  Figure  C5.  The  use  of  t*  =  0.1  sec 
for  Pg  may  not  be  justified,  but  a  different  value  will  simply  raise  or  lower  all  data  points  by 
the  same  amount  so  that  the  correlation  coefficient  and  mean  slope  values,  indicated  for  the 
three  plots  in  Figure  C5,  will  not  change.  The  linear  trends  (dashed  lines)  have  correlation  ! 

coefficients  of  0.846,  0.977,  and  0.992,  respectively,  again  suggesting  the  Pg  spectra  are  , 

I 

shaped  the  most  by  ESD,  somewhat  less  by  scaled  depth,  and  the  least  by  shot  depth.  A  com¬ 
parison  of  the  mean  slope  values  for  the  data  in  Figures  4  and  5  (indicated  for  each  plot)  indi¬ 
cates  Pn  spectral  slope  to  be  more  variable  than  the  Pg  spectra  slope.  i 

I 

i 

An  examination  of  the  results  in  Figures  1  through  5  shows  internal  consistency  and  the  I 

t 

1 

following  two  results:  (1)  spectral  ratios  Pn/Lg,  Pg/Lg  as  well  as  the  spectra  of  Pn  and  Pg 

i 

i 

l 
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Figure  C4.  Mean  spectral  slope  for  Pn  (corrected  for  noise,  yield,  and  t*)  over  the  frequen¬ 
cy  range  1  to  5  Hz  versus  (a)  log  (shot  depth  in  meters),  (b)  scaled  depth,  and  (c) 
log  ESD  for  7  explosions,  numbered  as  in  Table  Cl.  The  vertical  lines  represent 
error  bars  with  one  standard  deviation.  The  correlation  coefficient  and  mean  slope 
value  of  the  least  squares  linear  trend  (dashed  line)  is  indicated  for  each  plot. 
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depend  the  most  on  ESD,  somewhat  less  on  scaled  depth,  and  the  least  on  shot  depth,  and  (2) 
the  enrichment  of  higher  frequencies  due  to  increased  ESD  is  the  largest  in  the  spectra  of  Pn, 
somewhat  less  for  Pg,  and  the  least  for  Lg. 


Low-Frequency  Asymptote  and  Deconvolved  Pn 

The  low-frequency  asymptote,  A0,  and  the  area  under  the  initial  pulse  of  deconvolved  Pn 
(denoted  by  fi0)  are  good  measures  of  the  explosion  yield  (Douglas  et  al.,  1987;  McLaughlin 
et  al.,  1987).  A  spectral  measure  of  A0  was  made  by  removing  the  instrument  response, 
correcting  for  noise,  and  then  least-squares  fitting  the  log-amplitude  spectra  to  a  model  spec¬ 
trum  with  a  fall-off  rate  of  f-3  (McLaughlin  et  al.,  1987).  The  amplitude  spectrum  was 
assumed  to  have  the  form  (Boatwright,  1978) 


A(f)  = 


A0 

Vi  +(f/fc)6 


(C7) 


where  fc  is  the  comer  frequency.  The  least  squares  fit  to  the  log-amplitude  spectra  was  deter¬ 
mined  for  the  combination  of  A0  and  fc;  the  derived  values  are  listed  in  Table  C2.  Plots  of 
log  ( AJY )  versus  (a)  shot  depth,  (b)  scaled  depth,  and  (c)  ESD,  are  shown  in  Figure  C6. 
Assuming  linear  trends  (dashed  lines),  it  seems  that,  for  a  fixed  yield,  the  low-frequency 
asymptote,  A0,  is  influenced  the  most  by  scaled  depth,  somewhat  less  by  ESD,  and  consider¬ 
ably  less  by  shot  depth.  The  values  of  mean  slope  of  the  linear  trends  in  the  three  plots  in  Fig¬ 
ure  C6  are  also  indicated. 
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Figure  C6.  Log  (AJY)  versus  (a)  shot  depth,  (b)  scaled  depth,  and  (c)  log  ESD  for  7  ex¬ 
plosions,  numbered  as  in  Table  Cl.  The  correlation  coefficient  and  mean  slope 
value  of  the  least  squares  linear  trend  (dashed  line)  is  indicated  for  each  plot. 
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TABLE  C2 

VALUES  OF  A0,  fc,  AND  FOR  PAHUTE  MESA  EXPLOSIONS 


No. 

Name 

A0 

(nm-sec) 

fc 

(Hz) 

Oo 

(nm-sec) 

1 

BUTEO 

16.62 

6.64 

13 

2 

REX 

189.0 

3.15 

160 

3 

DURYEA 

310.0 

2.68 

414 

4 

SCOTCH 

936.9 

2.60 

1257 

5 

CABRIOLET 

12.77 

3.81 

11.5 

6 

SCHOONER 

103.0 

1.68 

142 

7 

BENHAM 

10760 

1.89 

8627 

The  Pn  arrivals  at  the  LRSM  station  KN-UT  were  deconvolved  by  removal  of  the  instru¬ 
ment  response  and  selection  of  an  optimum  bandpass  that  takes  into  consideration  the  ambient 
noise  spectrum  (Gupta  and  McLaughlin,  1987).  The  values  of  Qq  were  obtained  by  measuring 
the  area  of  the  first  positive  pulse  (Table  C2).  Figure  C7  shows  plots  of  log  (QJY)  versus  (a) 
shot  depth,  (b)  scaled  depth,  and  (c)  ESD.  The  three  plots  have  linear  trends  (dashed  lines) 
suggesting  that  £20  is  influenced  the  least  by  shot  depth  and  considerably  more  by  scaled  depth 
and  ESD.  The  values  of  mean  slope  of  the  linear  trends  in  the  three  plots  in  Figure  C7  are  also 
indicated;  these  are  all  smaller  than  the  corresponding  slopes  in  Figure  C6. 

As  expected,  the  two  measures  A0  and  Q0  are  nearly  equal  (Table  C2).  A  plot  (log-log 
scale)  of  the  area  under  the  initial  pulse  of  deconvolved  Pn,  Qq  versus  the  low-frequency 
asymptote,  A0,  is  shown  in  Figure  C8.  The  near-equivalence  of  the  two  measures  is  attested 
by  the  high  correlation  coefficient  and  an  estimated  slope  value  close  to  unity. 

When  comparing  contained  and  cratering  explosions,  Blandford  (1976)  observed 
significant  spectral  differences  between  the  two  types  of  explosions  for  both  body  and  surface 
waves.  It  is  therefore  important  to  examine  whether  scaled  depth  and  ESD  are  significant 
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Figure  Cl.  Similar  to  Figure  C6  for  log  (fl/V).  Note  that  the  mean  slope  values  of  the 
linear  trends  in  (a),  (b),  and  (c)  are  smaller  than  the  corresponding  values  in  Figure 
C6,  indicating  somewhat  smaller  variation  in  log  (£VY)  as  compared  to 
log  (Aq/Y). 


Figure  C8.  Plot  (log-log  scale)  of  the  area  under  the  initial  pulse  of  deconvolved  Pn,  £20 
versus  the  low-frequency  asymptote,  A0,  both  in  units  of  nm-sec.  The  neai 
equivalence  of  the  two  measures  is  confirmed  by  the  high  correlation  coefficien 
value  and  slope  close  to  unity. 
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parameters  when  the  two  cratering  shots  are  excluded  from  the  data  set.  The  data  in  Figures  C2 
and  C3  appear  to  suggest  poor  correlations  for  all  three  (shot  depth,  scaled  depth,  and  ESD) 
parameters  if  only  the  contained  shots  are  considered.  However,  the  data  in  the  remaining  four 
figures  (C4,  C5,  C6,  and  C7)  strongly  indicate  the  scaled  depth  and  ESD  to  provide  correla¬ 
tions  that  are  considerably  better  than  with  shot  depth.  A  comparison  of  the  correlation 
coefficient  values  for  all  7  shots  and  for  only  5  contained  shots  is  shown  in  Table  C3. 

TABLE  C3 

COMPARISON  OF  RESULTS  WITH  AND  WITHOUT  TWO  CRATERING  SHOTS 


CORRELATION  COEFFICIENT 

PARAMETER 

7  SHOTS 

5  SHOTS 

LOG  SHOT 

DEPTH 

LOG  SCALED 

DEPTH 

LOG  ESD 

LOG  SHOT 

DEPTH 

LOG  SCALED 

DEPTH 

LOG  ESD 

Pn  SPECTRAL  SLOPE 

0.827 

0.970 

0.979 

0.092 

0.908 

0.928 

Pg  SPECTRAL  SLOPE 

0.8^6 

0.977 

0.992 

0.145 

0.983 

0.987 

LOG(tVY) 

0.653 

0.978 

0.921 

0.133 

0.984 

0.994 

LOG  (AJY) 

0.545 

0.885 

0.819 

0.035 

0.883 

0.938 

Comparison  with  Mueller-Murphy  Scaling 

By  analyzing  a  large  amount  of  near -field  data,  Mueller  and  Murphy  (1971)  found  that 
the  yield,  depth  of  burial,  and  emplacement  medium  play  dominant  roles  in  defining  the 
seismic-spectrum  scaling  of  underground  nuclear  explosions.  They  used  source-receiver  dis¬ 
tances  of  about  100  km  or  less  and  examined  the  spectra  of  the  entire  seismograms.  In  nearly 
all  cases,  S  waves  (including  Lg)  dominate  over  other  arrivals.  In  fact,  a  comparison  of  the 
spectral  ratios  based  on  S  waves  and  total  seismograms  by  Murphy  (1975)  indicated  the  two 
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spectral  measures  to  be  essentially  identical  for  frequencies  greater  than  1  Hz. 

Our  results  clearly  indicate  that  the  spectral  scaling  is  strongly  dependent  on  the  phase 
type.  In  other  words,  the  variation  with  yield  of  the  source  spectrum  is  not  the  same  for  Pn, 
Pg,  or  Lg  phases.  For  example,  in  Figure  Cl,  the  inter-event  spectral  ratios  of  Pn,  Pg,  and  Lg 
are  significantly  different  for  both  pairs  of  shots.  In  most  cases,  S  waves  or  Lg  dominate  the 
ground  motion,  so  that  scaling  based  on  the  entire  seismogram  may  not  be  valid  for  scaling  of 
Pn  or  Pg.  As  an  example,  let  us  consider  the  data  from  Rex  (overburied)  and  Schooner 
(cratering),  both  of  which  were  used  in  Mueller  and  Murphy’s  study.  One  of  their  results  is 
reproduced  in  Figure  C9,  which  shows  the  pseudo  relative  velocity  (PSRV)  spectra  (for  5  % 
damping)  of  Rex  and  Schooner  based  on  the  horizontal  component  data  at  station  SE-6  (epi- 
central  distance  about  110  km).  We  obtained  the  PSRV  response  spectrum  in  the  manner 
described  in  Gupta  and  McLaughlin  (1987);  only  spectral  points  with  S/N  power  ratio  of  at 
least  2  are  plotted.  The  PSRV  spectra  from  the  horizontal  radial  (R)  component  data  at  KN-UT 
are  shown  in  Figure  CIO,  which  shows  results  from  both  (a)  Pn  (6.4  sec  window)  and  (b)  total 
seismogram  (denoted  by  ALL).  There  is  good  agreement  between  the  spectra  in  Figures  9  and 
10b,  especially  for  periods  longer  than  0.2  sec.  In  both  figures,  the  Schooner  and  Rex  spectra 
cross  each  other  at  a  period  of  about  0.5  sec.  The  PSRV  spectra  of  Pn  (Figure  ClOa)  appear  to 
be  quite  different  from  those  in  Figure  Cl  Ob,  especially  for  periods  longer  than  0.5  sec.  The 
spectral  ratios  Rex/Schooner  are  probably  not  significantly  different,  except  that  for  periods 
longer  than  0.5  sec,  the  ratio  is  about  2  for  Pn  and  about  3  for  ALL  in  Figure  CIO.  The  results 
from  another  pair  of  events,  Buteo  (overburied)  and  Cabriolet  (cratering),  are  shown  in  Figure 
Cl  1.  The  spectra  based  on  the  total  seismogram  (Figures  Cl  Ob  and  Cl  lb)  appear  fairly  similar, 
whereas  those  for  Pn  (Figures  ClOa  and  Cl  la)  are  significantly  different  For  periods  greater 
than  about  0.5  sec,  the  spectral  ratio  Buteo/Cabriolet  for  Pn  (Figure  Cl  la)  is  significantly 
different  from  that  for  the  total  seismogram  (Figure  Cl  lb). 
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Figure  Cll.  Similar  to  Figure  CIO  for  Buteo  (overburied)  and  Cabriolet  (underburied). 
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Quite  often,  horizontal  component  data  are  not  available  and  only  vertical  component 
records  are  used.  On  the  vertical  component  record  of  Rex,  the  phase  Pg  is  clipped,  so  that 
PSRV  spectra  similar  to  those  in  Figure  ClOb  could  not  be  obtained.  Results  from  the  vertical 
component  data  from  Buteo  and  Cabriolet  at  KN-UT  are  shown  in  Figure  Cl 2.  A  comparison 
with  Figure  C 1 1  suggests  that  the  horizontal  and  vertical  component  spectra  are  by  no  means 
identical,  although  the  spectral  ratios  Buteo/Cabriolet  are  rather  similar. 

Discussion  and  Conclusion 

The  spectral  ratio  P/P-coda,  obtained  from  short-period,  vertical  component  records  of  20 
NTS  explosions  recorded  at  teleseismic  distances,  showed  significant  correlation  with  both  P- 
wave  overburden  velocity,  a,  and  shot  depth  (see  Figures  4  and  5,  Gupta  and  Blandford, 
1987).  Assuming  that  these  shots  were  fired  approximately  at  their  scaled  depths,  equation  (2) 
would  suggest  that  both  a  and  h  are  rough  measures  of  ESD,  since  a  and  p  are  nearly  always 
well  correlated  and  each  increases  with  depth.  Thus,  the  spectral  slope  of  P/P-coda  at  telese¬ 
ismic  distances  varies  with  ESD  in  the  same  manner  as  Pn/Lg  at  regional  distances.  Moreover, 
in  Gupta  and  Blandford’s  (1987)  study,  most  of  the  spectral  variation  appeared  to  be  in  the  ini¬ 
tial  P  waves  and  P-coda  was  found  to  be  stable,  similar  to  the  results  for  Pn  and  Lg  spectra  in 
this  study.  By  applying  frequency-wavenumber  methods  to  NORESS  data  from  NTS  explo¬ 
sions,  Dainty  (1985)  found  evidence  indicating  that  the  dominant  mode  of  near-source  scatter¬ 
ing  contributing  to  teleseismic  P  coda  is  Lg  to  P  due  to  shallow  heterogeneity.  It  seems  there¬ 
fore  that  the  spectral  variations  with  ESD  in  Pn  and  Lg  at  regional  distances  are  similar  to 
those  of  P  and  P-coda  at  teleseismic  distances. 

Flynn  (1986)  studied  effects  of  source  depth  on  near  source  seismograms  by  examining 
the  spectra  of  P  and  SV-Rayleigh  surface  waves  from  five  115  kg  chemical  explosions  buried 
at  depths  ranging  from  1.8m  (optimum  cratering  depth)  to  11.5  m  (fully  contained).  Energy  in 
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Figure  C12.  Similar  to  Figure  CIO  for  Buteo  (overburied)  and  Cabriolet  (underburied)  but 
for  the  vertical  (Z)  components. 
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the  high  frequency  band  of  10  to  40  Hz,  dominated  by  P  waves,  was  found  to  increase 
significantly  with  source  depth  whereas  energy  in  the  low  frequency  band  of  4  to  12  Hz,  con¬ 
sisting  mostly  of  SV-Rayleigh  waves,  was  not  a  strong  function  of  shot  depth.  It  was  con¬ 
cluded  that  the  fully  contained  source  generates  more  high  frequency  energy  as  a  result  of  an 
increase  of  coupling  over  the  near-surface  source.  Our  results  regarding  the  influence  of  scaled 
depth  (or  ESD)  on  the  relative  excitation  of  P  and  Lg  as  well  as  on  A0  and  Qq  are  therefore  in 
good  agreement  with  those  of  Flynn  (1986). 

Underground  nuclear  explosions  are  generally  very  complex  with  both  primary  and 
secondary  sources  of  seismic  waves  (e.g.  Stump,  1987;  Stump  and  Reinke,  1987).  The  ideali¬ 
zation  of  an  explosion  source  as  a  pressure  pulse  applied  to  the  walls  of  a  spherical  cavity 
should  at  best  be  used  to  infer  only  qualitative  results.  Non-linear  effects,  known  to  be  impor¬ 
tant  for  most  nuclear  explosions  (especially  those  with  shallow  scaled  depths)  also  drastically 
modify  the  seismic  signals  (Perl  et  al.,  1979;  Rimer  et  al.,  1979;  Day  et  al.,  1986).  Note  that 
non-linear  finite  difference  modeling  of  the  explosion  source  indicates  strong  dependence  on 
scaled  depth  that  differs  significantly  from  that  derived  by  assuming  an  elastic  point  source 
(Day  et  al.,  1986). 

Pn  arrivals  from  deeper  and/or  overburied  explosions  are  richer  in  high  frequencies  as 
compared  to  their  Lg.  These  results,  similar  to  the  earlier  results  of  Blandford  (1976)  based  on 
observations  of  Pn  and  Rayleigh  wave  speed  a,  indicate  that  the  scaling  laws  may  be  strongly 
dependent  on  the  phase  type.  Scaled  depth  or  ESD,  a  measure  of  the  source  size,  appears  to 
play  an  important  role  in  shaping  the  spectra  of  various  regional  phases  as  well  as  in  determin¬ 
ing  A0  and  ft0,  both  of  which  are  measures  of  the  explosion  seismic  moment 


Teledyne  Geotech 


52 


November  1987 


ANNUAL  REPORT 


STUDIES  IN  DECOUPLING 


TGAL-87-4 


D.  DECONVOLUTION  AND  SPECTRAL  RATIOS  FROM  SOVIET  SALT  SHOTS 

Introduction 

We  first  applied  the  maximum-likelihood  multichannel  deconvolution  method  of  Shum- 
way  and  Der  (1985)  to  teleseismic  P  waves  from  10  salt  shots  in  the  Astrakhan  region  of  the 
USSR  (Table  Dl).  The  deconvolution  method  utilizes  the  well  known  fact  that  P  wave  spectra 
can  be  decomposed  into  source  and  receiver  spectral  factors;  details  of  the  procedure  can  be 
found  in  Shumway  and  Der  (1985)  and  Der  et  al.  (1987).  The  multichannel  deconvolution 
approach  for  estimating  source  and  site  region  characteristics  has  several  advantages.  First,  the 
method  utilizes  both  amplitude  and  phase  information  in  the  frequency  domain  calculations. 
Since  this  is  a  frequency  domain  calculation,  the  maximum  available  bandwidth  can  be  used, 
avoiding  the  ambiguities  associated  with  time  domain  fitting  of  narrow  frequency  band  seismic 
signals.  Furthermore,  no  a  priori  assumptions  need  to  be  made  about  the  complexity  of  either 
the  source  or  the  site  spike  series,  e.g.  the  presence  of  pP  does  not  need  to  be  presupposed  by 
providing  an  initial  guess  of  the  pP  delay  time  and  amplitude.  Therefore,  the  source  functions 
obtained  in  the  deconvolution  process  provide  a  better  estimate  of  the  explosion  source  time 
functions,  especially  the  presence  or  lack  of  secondary  arrivals  following  the  direct  P  wave 
such  as  pP  or  spall. 

We  also  carried  out  single-station  deconvolutions  for  4  salt  shots  of  October  16,  1982 
recorded  at  two  SRO  stations,  ANMO  and  KONO.  In  order  to  explore  why  the  last  shot  in  the 
sequence  has  a  larger  magnitude  value  than  the  earlier  three,  we  obtained  mean  amplitude 
ratios  of  P  and  P-coda,  spectral  ratio  slopes  of  P,  and  P/P-coda  spectral  slopes. 
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Multichannel  Deconvolutions 


Two  sequences  of  Astrakhan  events  containing  4  and  6  events  on  October  16,  1982  and 
September  24,  1983,  respectively  (Table  Dl),  were  deconvolved  from  observations  at  GBA  and 
the  four  RSTN  stations,  RSON,  RSON,  RSNT,  and  RSSD.  From  the  GBA  array  data,  four 
common  receivers  with  good  S/N  were  selected.  The  values  of  t*  for  GBA  and  RSTN  were 
assumed  to  be  0.15  and  0.2  sec,  respectively.  Following  Borg  (1983),  the  yields  of  all  10  shots 
in  salt  were  considered  to  be  15  kt  each;  the  corresponding  von  Seggem  and  Blandford  (1972) 
or  "VSB"  wavelet  was  based  on  this  yield  value. 

TABLE  Dl 

ASTRAKHAN  SHOTS  DECONVOLVED  IN  THE  STUDY 
(Information  from  the  ISC  Bulletins) 


Event 

Origin  Time 

Location 

'"b 

Lat(°N) 

Lon(°E) 

1982289a 

05:59:57.4 

46.77 

48.22 

5.2 

1982289b 

06:04:57.4 

46.77 

48.24 

5.2 

1982289c 

06:09:57.4 

46.77 

48.22 

5.2 

1982289d 

06:14:57.5 

46.75 

48.20 

5.4 

1983267a 

04:59:57.2 

46.81 

48.29 

5.2 

1983267b 

05:04:57.3 

46.81 

48.27 

5.1 

1983267c 

05:09:57.4 

46.84 

48.33 

5.0 

1983267d 

05:14:57.4 

46.84 

48.29 

5.2 

1983267e 

05:19:57.4 

46.88 

48.30 

5.4 

1983267f 

05:24:57.0 

46.80 

48.29 

5.3 

Figure  Dl  shows  source  time  functions  for  the  four  events  of  October  16,  1982  decon¬ 
volved  separately  from  data  recorded  at  GBA  (left)  and  four  of  the  RSTN  stations  (center). 
The  VSB  wavelet  has  been  removed  in  the  deconvolutions.  The  relatively  great  depth  of  these 
shots  can  be  inferred  from  the  pP  delay  time  of  over  1  sec  for  all  four  shots.  There  are  also 
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two  or  three  arrivals  which  consistently  show  up  between  the  P  and  pP.  These  appear  to  be 
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similar  to  those  observed  for  Salmon,  and  are  interpreted  to  be  reflections  from  the  high 
impedance  contrast  layers  between  the  shot  point  and  the  surface  (Blandford  and  McElfresh, 
1982).  This  is  a  distinct  possibility  because  these  shots  were  detonated  within  an  area  of  salt 
domes. 

The  source  time  functions  on  the  right  in  Figure  D1  are  the  result  of  deconvolving  the 
two  separate  sets  of  deconvolutions  at  the  left  and  center.  In  this  case,  t*  and  instrument 
response  were  not  removed  since  they  had  already  been  removed  in  the  previous  deconvolu¬ 
tions.  These  "joint"  deconvolutions  contain  parts  of  the  source  time  functions  that  are  common 
to  observations  at  both  arrays  and  are  therefore  somewhat  more  azimuthally  averaged  and  less 
complex  than  the  other  two  sets.  However,  they  do  still  show  clear  P  and  pP  phases,  which 
incidentally  look  very  much  like  mirror  images  of  each  other. 

Figure  D2  shows  the  deconvolved  source  time  functions  for  the  six  events  of  September 
24,  1983  as  recorded  at  GBA  and  obtained  with  the  same  parameters  as  for  the  four  shots  in 
Figure  Dl.  As  with  the  previous  set  of  events,  the  pP  delay  time  is  over  1  sec,  and  there  is  an 
additional  arrival  or  two  between  the  P  and  pP.  The  six  events  are  very  closely  co-located  in 
space  and  time,  so  it  is  probably  not  surprising  that  their  source  functions  look  so  similar. 

Figure  D3  shows  the  deconvolved  source  time  functions  for  all  ten  Astrakhan  events,  as 
recorded  at  GBA,  deconvolved  together.  Results  from  the  RSTN  stations,  with  the  same 
parameters  as  for  the  four  shots  in  Figure  Dl,  are  shown  in  Figures  D4  and  D5.  One  shot  was 
not  recorded  at  the  RSTN  station  RSSD.  Source  deconvolutions  from  all  ten  shots  recorded  at 
three  RSTN  stations  (Figure  D4)  are  similar  to  those  for  9  shots  based  on  data  from  the  four 
RSTN  stations  (Figure  D5).  Again,  the  depth  and  other  secondary  arrivals  are  very  similar  for 
all  events. 
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Figure  D2.  Deconvolved  source  time  functions  for  the  six  Astrakhan  events  of  September 
24,  1983  recorded  at  GBA.  The  arrows  denote  the  peaks  of  P  and  pP  arrivals. 


1982289d 


4.0  SEC 


Figure  D3.  Deconvolved  source  time  functions  for  the  ten  Astrakhan  events  of  October  16, 
1982  and  September  24,  1983  recorded  at  GBA. 
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Figure  D4.  Deconvolved  source  time  functions  for  the  ten  Astrakhan  events  of  October  16, 
1982  and  September  24,  1983  recorded  at  the  three  RSTN  stations  RSON,  RSNT, 
and  RSNY. 
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Figure  D5.  Deconvolved  source  time  functions  for  nine  Astrakhan  events  (four  on  October 
16,  1982  and  five  on  September  24,  1983)  recorded  at  the  four  RSTN  stations 
RSON,  RSNT,  RSNY,  and  RSSD. 
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Single  Station  Deconvolutions 

Deconvolution  of  the  four  Astrakhan  explosions  of  October  16,  1982  was  extended  to  the 
teleseismic  data  from  the  SRO  stations  ANMO  and  KONO  at  epicentral  distances,  A  of 
approximately  96°  and  26°,  respectively.  The  deconvolution  is  accomplished  by  removal  of  the 
instrument  response,  correction  for  the  appropriate  t*,  and  selection  of  an  optimum  bandpass 
that  takes  into  consideration  both  the  instrument  response  and  the  ambient  noise  spectra  (Gupta 
and  McLaughlin,  1987).  We  used  t*  =  0.45  for  ANMO  and  t*  =  0.40  for  KONO.  The  results 
are  shown  in  Figure  D6.  The  first  P  arrivals  at  KONO  are  weak  but,  because  of  the  triplication 
in  the  travel-time  curves  for  the  epicentral  distance  to  this  station,  a  second  set  of  P  waves 
corresponding  to  a  later  branch  of  the  travel-time  curve  can  be  seen.  The  results  are  similar  to 
those  obtained  earlier  from  the  RSTN  stations  and  the  GBA  array;  the  pP  and  P  arrivals  are 
again  separated  by  about  1  sec  for  each  of  the  four  shots. 

Spectral  Ratios  of  P  Waves 

Deconvolved  waveforms  of  the  four  shots  of  October  16,  1982  are  similar  and  the  P  and 
pP  arrivals  are  separated  by  nearly  the  same  amount,  suggesting  nearly  the  same  depth  for  all 
shots.  Yet  the  magnitude  of  the  last  shot  in  the  sequence  is  0.2  unit  larger  than  those  of  the 
earlier  shots  (Table  Dl).  In  order  to  explore  a  possible  reason  for  this,  we  examined  the  spectra 
of  P  and  P-coda  for  all  four  shots.  We  first  obtained  the  mean  spectral  amplitude  ratios  of 
6.4-sec  long  P  windows  (with  10%  cosine  taper  and  corrected  for  noise),  averaged  over  the 
frequency  range  of  1  to  5  Hz  (the  frequency  range  with  generally  good  S/N),  for  all  inter-event 
combinations.  Results  for  the  six  possible  inter-event  combinations  from  4  shots  (denoted  by 
El,  E2,  E3,  and  E4,  respectively)  recorded  at  the  four  RSTN  stations  and  four  sensors  of  the 
GBA  array  are  shown  in  Figure  D7.  On  the  average,  the  last  shot  in  the  sequence  is  about  0.2 
magnitude  unit  larger  than  the  first  three,  in  agreement  with  their  mb  values.  Similar  mean 


Teledyne  Geotech 


61 


November  1987 


2»  0  SEC  12H1.0  Nti 


2.0  SEC  1931.5  NM  ASTRAKHAN  AT  KONO  (A  =  26°) 

Figure  D6.  Deconvolved  source  time  functions  for  the  four  Astrakhan  events  of  October 
16,  1982  recorded  at  ANMO  and  KONO.  The  numbers  denote  zero-to-peak  ampli¬ 
tudes  in  nanometers.  The  arrows  denote  the  peaks  of  P  and  pP  arrivals. 
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Figure  D7.  Inter-event  mean  spectral  amplitude  ratios  of  P,  in  log  units  and  over  the  fre¬ 
quency  range  1  to  5  Hz,  for  4  RSTN  and  4  GBA  sensors.  The  ratios  indicate  the 
fourth  event  to  be,  on  the  average,  larger  than  the  first  three  by  about  0.2  unit. 
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amplitude  ratios  derived  from  the  spectra  of  12.8  sec  of  P  coda  (again  with  10%  cosine  taper 
and  corrected  for  noise),  immediately  following  the  6.4  sec  long  P  window,  are  shown  in  Fig¬ 
ure  D8.  The  P  coda  of  the  last  shot  appears,  on  the  average,  to  be  about  0.1  magnitude  unit 
larger  than  those  of  the  first  three.  A  comparison  of  Figures  D7  and  D8  shows  similar  trends 
but  generally  smaller  variation  for  P  coda  than  for  P. 

Inter-event  mean  spectral  slopes  of  P,  with  10%  cosine  taper  and  corrected  for  noise,  over 
the  frequency  range  of  1  to  5  Hz,  at  both  RSTN  and  GBA,  are  shown  in  Figure  D9.  Each  of 
the  spectral  ratios  E4/E1,  E4/E2,  and  E4/E3  has  a  mean  value  of  about  -0.1/Hz.  This  "red- 
shift"  in  the  spectrum,  of  the  last  shot  is  significantly  larger  than  that  expected  on  the  basis  of 
the  0.2  magnitude  unit  difference  in  magnitudes  (the  latter  value  is  about  -0.02/Hz  on  the  basis 
of  theoretical  spectra  with  yields  estimated  from  the  mb  values).  Thus  the  last  event  appears  to 
be  rich  in  low  frequency  (or  depleted  in  high  frequencies)  as  compared  to  the  first  three  shots. 

We  also  obtained  the  P/P-coda  spectral  slopes  for  the  four  shots  at  the  four  RSTN  sta¬ 
tions.  Results  for  the  mean  slopes,  over  the  frequency  range  of  1  to  5  Hz,  are  shown  in  Figure 
DIO.  At  each  station,  the  spectral  slope  value  for  the  last  shot  is  lower  than  for  the  first  three 
shots.  Results  from  the  GBA  array  were  found  to  be  unreliable  due  to  poor  S/N. 

In  a  comparison  of  P  spectra  from  several  shots  in  close  proximity  to  one  another,  Aki  et 
al.  (1970)  observed  the  later  shots  to  have  larger  amplitudes  and  lower  frequency  content  than 
the  first  shot.  They  attributed  this  to  changes  in  the  medium  properties  due  to  shock  waves 
from  the  earlier  shot;  the  inelastic  zone  from  the  first  shot  contained  crushed  or  cracked  rocks 
which,  for  a  later  shot,  preferentially  attenuated  the  high  frequencies  and  also  gave  rise  to 
greater  displacement  Our  results  in  Figures  D7  and  D9  are  therefore  consistent  with  the  possi¬ 
bility  that  the  fourth  shot  was  detonated  within  the  inelastic  zone  of  an  earlier  shot.  Gupta  and 
Blandford  (1987)  compared  the  spectra  of  P  and  P-coda  at  teleseismic  distances  for  shots  in 
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Figure  D9.  Inter-event  mean  spectral  ratio  slopes  of  P,  in  log  units  per  Hz  and  over  the 
frequency  range  1  to  5  Hz,  for  4  RSTN  and  4  GBA  sensors.  The  slope  values  in¬ 
dicate  the  fourth  event  to  have,  on  the  average,  about  O.i  unit/Hz  lower  slope  than 
the  first  three  shots. 
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Figure  DIO.  Mean  P/P-coda  spectral  slopes,  in  log  units  per  Hz  and  over  the  frequency 
range  1  to  5  Hz,  for  the  four  shots  of  October  16,  1982  recorded  at  the  four  RSTN 
stations.  The  vertical  lines  represent  error  bars  with  one  standard  deviation.  At 
each  station,  the  fourth  shot  has  a  lower  spectral  slope  value  than  for  the  first  three 
shots. 
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various  velocity  media  and  found  both  amplitude  and  spectral  variation  in  P  coda  to  be  consid¬ 
erably  smaller  than  in  P,  similar  to  the  results  in  Figures  D7  and  D8.  Furthermore,  they  found 
lower  P/P-coda  spectral  slope,  as  for  the  fourth  shot  in  Figure  DIO,  to  be  an  indication  of 
lower  medium  velocity  around  the  shot  point.  It  is  therefore  quite  possible  that  the  yield  of  the 
fourth  snot  was  no  larger  than  the  three  earlier  shots  but  that  it  was  in  close  proximity  of  one 
or  more  earlier  shots. 

According  to  Borg  (1983),  the  four  explosions  occurred  within  five  minutes  of  one  other, 
possibly  to  minimize  disturbance  at  the  city  of  Astrakhan  some  30  to  40  km  distant,  and  the 
yields  of  these  shots  were  about  15  kt  each.  The  Soviets  are  known  to  detonate  shots  near  the 
cavity  of  a  previous  shot  (e.g.  Nordyke,  1973,  Figure  30;  Polskov  et  al.,  1980,  Figure  2). 
Therefore,  it  seems  likely  that  the  four  shots  were  of  equal  yields,  fired  close  to  each  other, 
and  the  last  shot  was  influenced  by  the  fracture  zone  of  an  earlier  shot 


Teledyne  Geotech 


68 


November  1987 


ANNUAL  REPORT 


STUDIES  IN  DECOUPLING 


TGAL-87-4 


E.  ANALYSIS  OF  DATA  FROM  SOVIET  SALT  SHOTS  WITH  KNOWN  YIELDS 

We  analyzed  teleseismic  P  arrivals  from  three  USSR  salt  explosions  with  known  yields 
and  shot  depths  (Blandford,  1978)  as  recorded  at  the  LRSM  stations  NP-NT  and  RK-ON.  Two 
of  these  (shots  of  22  April  1966  and  1  July  1968)  had  their  epicenters  in  the  Azgir  region, 
approximately  48°N,  48°E,  north  of  the  Caspian  region.  The  epicenter  of  the  third  (21  May 
1968)  was  at  Bukara,  approximately  39°N,  65°E.  All  three  explosions  were  well  recorded  at 
the  three-component  stations  NP-NT  and  RK-ON.  The  vertical  component  records  of  these 
three  salt  shots,  along  with  the  U.S.  salt  shot  Salmon,  are  shown  in  Figures  El  and  E2.  The 
analog  record  of  the  1  July  1968  shot  at  both  RK-ON  and  NP-NT  was  redigitized,  as  the  ear¬ 
lier  waveforms  suffered  from  spikes  acquired  during  the  analog  to  digital  conversion. 

An  examination  of  the  spectra  of  the  three  Soviet  salt  explosions  with  known  yields  was 
made  in  order  to  test  whether  they  are  consistent  with  conventional  scaling  for  salt,  such  as 
von  Seggern  and  Blandford’s  (1972)  scaling  (equations  C3,  C4,  and  C5).  Of  the  three  Soviet 
shots,  only  two  (22  April  1966  and  1  July  1968)  were  very  near  each  other  so  that  the  propa¬ 
gation  path  effects  may  be  considered  to  be  the  same.  The  spectral  ratios  (I  Jul  68)/(22  Apr 
66)  for  P  (6.4  sec  windows)  and  P  coda  (12.8  sec  windows)  at  both  NP-NT  and  NP-NT  are 
shown  in  Figure  E3.  The  spectra  were  obtained  with  10  %  cosine  taper  and  no  smoothing.  The 
spectral  ratios  were  corrected  for  noise  and  only  points  with  S/N  power  ratio  of  at  least  2  are 
shown.  The  theoretical  spectral  ratios,  also  shown  in  Figure  E3,  were  obtained  by  using  values 
of  the  two  medium-dependent  constants,  kg  and  B  to  be  26  and  0,  respectively. 

Lyuke  et  al.  (1976)  used  the  von  Seggern  and  Blandford  source  model  to  scale  success¬ 
fully  the  spectra  of  salt  shots  recorded  at  several  distances.  Their  analysis,  which  included  the 
1.1  kt  shot  of  22  April  1966,  indicated  the  best  agreement  between  theory  and  observation 
when  k0  and  B  were  selected  to  be  11  and  0,  respectively.  They  also  considered  the  values  lq, 


Teledyne  Geotech 


69 


November  1987 


4  SALT  SHOTS  RK-ON 

Figure  E2.  Records  of  the  same  four  salt  shots  as  in  Figure  El  at  the  LRSM  station  RK- 
ON.  A  denotes  the  epicentral  distance  in  degrees. 
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=  11,  B  =  1.  Theoretical  spectral  ratios  for  these  two  sets  of  values  of  k,,  and  B  are  also  shown 
in  Figure  E3.  Considering  all  four  spectral  ratios,  the  theoretical  spectra  with  ko  =11  and  B  = 
0  appear  to  be  the  closest  to  the  observed  data.  It  seems  therefore  that  for  the  Soviet  salt  shots, 
k0  is  considerably  smaller  than  for  the  U.S.  salt  shot  Salmon.  Since  the  explosion  rise  time 
varies  inversely  as  k*,,  the  Soviet  salt  is  indicated  to  be  "softer"  than  the  Salmon  salt.  This 
could  be  due  to  the  Salmon  salt  being  harder  (possibly  because  of  larger  depth  and  overburden 
pressure)  than  the  USSR  salt  shot  as  suggested  by  a  comparison  of  their  physical  properties, 
listed  in  Table  El. 

TABLE  El 

COMPARISON  OF  SALMON  AND  USSR  1.1  KT  SHOT  SALT  PROPERTIES 


Property 

SALMON 

USSR  Shot 

Yield 

Depth 

Scaled  depth 
P-wave  velocity 
Density 

Young’s  modulus 
Reference 

5.3  kt 

828  m 

475  m/kt1/3 
4.67  km/sec 

2.2  gm/cc 

4.0  x  105  kb 
Mueller  (1969) 

1.1  kt 

161  m 

156  m/kt1/3 

4.2  km/sec 

2.15  gm/cc 

3.2  x  105  kb 
Lyuke  et  al.  (1976) 

We  used  the  maximum-likelihood  multichannel  deconvolution  method  of  Shumway  and 
Der  (1985),  described  in  Section  D,  to  obtain  the  deconvolved  source  terms  for  the  three  USSR 
salt  explosions  with  known  yields.  Figure  E4  shows  the  deconvolved  waveforms  which  are 
obtained  by  making  simultaneous  use  of  the  vertical  and  radial  component  records  at  the 
LRSM  station  NP-NT.  The  source  functions  have  been  corrected  for  instrumental  response  and 
attenuation  (by  assuming  t*  =  0.2),  but  no  correction  due  to  differences  in  yields  has  been 
applied.  Note  that  the  scaled  depths  (denoted  by  S.D.  in  Figure  E4)  vary  considerably.  Only 
the  two  deeper  (and  larger  scaled  depth)  explosions  show  clear  evidence  of  P  and  pP  arrivals; 
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Figure  E4.  Deconvolved  source  functions  of  three  USSR  salt  shots  derived  from  their  vert¬ 
ical  and  radial  component  records  at  NP-NT.  The  arrows  denote  the  peaks  on  P 
and  pP  arrivals.  The  letters  W,  h,  S.D.,  and  P/P-coda  denote  yield  (kt),  shot  depth, 
scaled  depth,  and  average  value  over  the  frequency  range  of  0.5  to  3.0  Hz,  respec¬ 
tively. 
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the  pP  -  P  times  agree  well  with  the  expected  values  (Blandford,  1978).  The  average  ratio 
P/P-coda  (over  the  frequency  range  of  0.5  to  3  Hz)  values  are  also  indicated,  and  they  seem  to 
be  larger  for  larger  scaled  depths,  perhaps  because  an  overburied  shot  is  likely  to  generate  less 
near-source  coda  than  an  underburied  shot.  Note  that  reciprocal  of  the  average  ratio  P/P-coda 
is  a  measure  of  the  complexity  (e.g.  Aki,  1982). 
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